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"...the most beautiful and important 
discovery of Mechanics." 

Lagrange to Maupertuis (November 1756) 



Abstract 

We review the development and practical uses of a generalized 
Maupertuis least action principle in classical mechanics, in which the 
action is varied under the constraint of hxed mean energy for the trial 
trajectory. The original Maupertuis (Euler-Lagrange) principle con- 
strains the energy at every point along the trajectory. The generalized 
Maupertuis principle is equivalent to Hamilton's principle. Recipro- 
cal principles are also derived for both the generalized Maupertuis and 
the Hamilton principles. The Reciprocal Maupertuis Principle is the 
classical limit of Schrodinger's variational principle of wave mechan- 
ics, and is also very useful to solve practical problems in both classical 
and semiclassical mechanics, in complete analogy with the quantum 
Rayleigh-Ritz method. Classical, semiclassical and quantum varia- 
tional calculations are carried out for a number of systems, and the 
results are compared. Pedagogical as well as research problems are 
used as examples, which include nonconservative as well as relativistic 
systems. 
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1 Introduction and History 



Variational Principles have a long and distinguished history in Physics. Apart 
from global formulations of physical principles, equivalent to local differential 
equations, they also are useful to approximate problems too difficult for an- 
alytic solutions. In the past century quantum variational calculations have 
been ubiquitous as an approximate method for the ground state of many 
difficult systems pQ. In recent years there has been some progress in using 
classical variational principles (Action Principles) to approximate the motion 
of classical systems such as classical molecules. We review recent develop- 
ments in this area of Classical Mechanics, although we also overlap with 
Quantum Mechanics. In particular we shall also discuss the use of the quan- 
tum variational principle for excited states, and the connection to classical 
action principles. 

In Classical Mechanics, variational principles are often called Least Ac- 
tion Principles, because the quantity subject to variations is traditionally 
the Action. To confuse matters there are two classical Actions, correspond- 
ing to two main Action Principles, which are called respectively, Hamilton's 
Action (S) and Maupertuis' Action (W); these differ from each other (they 
are related by a Legendre transformation), and the notation S, W is not 
universal (we follow current usage but in fact some authors use the same 
symbols reversed [2]). Both Actions, S and W, have the same dimensions 
(i.e. energyxtime, or angular momentum). Maupertuis' Least Action Princi- 
ple is the older of the two (1744) by about a century, and as we shall discuss, 
it is inconveniently formulated in all textbooks [3]. For clear statements of 
these old action principles, see the textbooks by Arnold [3], Goldstein et al 

and Sommerfeld We shall review these Action Principles in the next 
section. Some mathematics and notation required to deal with variational 
problems are summarized in Appendix I. The distinguishing feature of varia- 
tional problems is their global character j7j. One is searching for the function 
which gives a minimum (or stationary value) to an integral, as opposed to 
local extrema, where one is searching for the value of the variable which min- 
imizes ( or makes stationary) a function. Of course this dichotomy is not so 
sharp in practice, as the global problem is equivalent to a differential equa- 
tion (the Euler-Lagrange equation) which is of course local. The relation has 
even found its way into modern parlance "Think globally, act locally!". 

In Quantum Mechanics, where (usually) the aim is for approximations to 
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the energy by choosing a wave-function, one speaks of trial wave-functions, 
which are optimized according to a Variational Principle. In Classical Action 
Principles, the equivalent notion is that of a virtual path or trajectory, which 
should be optimized in accordance with the Action Principle. For uniformity 
we shall talk of Trial Trajectories in the classical domain, to emphasize the 
similarity of these ideas. 

As noted above, Maupertuis' Action Principle was formulated in precise 
form first by Euler and Lagrange some two and a half centuries ago; for 
the history, see the books by Yourgrau and Mandelstam and by Terral 0. 
However, even the formulation of these eminent mathematicians remained 
problematic, as emphasized by Jacobi, who stated in his Lectures on Dy- 
namics: " In almost all text-books, even the best, this Principle is presented 
so that it is impossible to understand" (TU1 HI] . To understand at least in 
part the problem (see [T5] for a detailed discussion), it is sufficient to note 
that Maupertuis' Principle assumes conservation of energy, whereas a well 
formulated principle, like Hamilton's Principle, implies energy conservation. 
As we shall see, it is possible to reformulate Maupertuis' Principle in a more 
general form to remedy this problem and also make it more useful in appli- 
cations OHl]- It should be noted that Hamilton's Principle does not suffer 
from any such drawbacks and is useful both conceptually (to derive equations 
of motion) and as a tool for approximations. The reformulated Maupertuis 
Principle has the advantage of being closely related to the classical limit 
of Schrodinger's Variational Principle of wave mechanics, and thereby lends 
itself easily to semiclassical applications. 

Even though we emphasize two Action Principles for Classical Mechanics, 
it is worth noting that there are many other formulations and reformulations 
of these principles. For a list of references on this point see ref.^2]- We note in 
particular the elegant work of Percival in the 1970's on variational principles 
for invariant tori (see ref. for a review). The emphasis on Maupertuis and 
Hamilton reflects the two main approaches of mechanics, one based on the 
Lagrangian L and the other based on the Hamiltonian H . 

In the above and in section 2 we use the traditional terminology for the 
variational principles, i.e. "least" action principles, but, since Jacobi's work 
it has been recognized jlj that the action is in fact stationary in general 
for the true trajectories. This means the first-order variation vanishes, and 
the action may be a minimum, a maximum or a saddle, depending on the 
second-order variation Hence "stationary" action principles would be 
more accurate terminology. Similar remarks apply to the general Maupertuis 
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and reciprocal variational principles discussed in section 2. 

In sections 2-4 we restrict ourselves to conservative holonomic systems. 
Various problems are solved using classical and semiclassical variational meth- 
ods and some comparisons are made with the results from the quantum varia- 
tional method for excited states ^E]- I n section 5 we discuss nonconservative 
(but nondissipative) systems, and in section 6 we mention briefly particular 
nonholonomic systems (which have velocity constraints), in connection with 
relativistic systems. The relations between classical and quantum variational 
principles (VP's) are discussed in section 7. 

2 Action Principles of Classical Mechanics 
2.1 Statements of Action Principles 

Hamilton's Least Action Principle (HP) states that, for a true trajectory of 
a system, Hamilton's Action S is stationary for trajectories which run from 
the fixed initial space-time point A = (g^t^) to the fixed final space-time 
point B = (q B ,t B ), 



where T = t B — tA is the duration. Here the (Hamilton) action S is the time 
integral of the Lagrangian L(q, q) from the initial point to the final point on 
the trial trajectory q(t), 



where q(t) is the generalized coordinate, q(t) = dq/dt the generalized velocity 
and t the time. In practice we choose tA = and t B = T for convenience. 
In general q stands for the complete set of independent generalized coor- 
dinates qi,q2,---,qf, where / is the number of degrees of freedom 17J. In 
(2.1) the constraint of fixed T is indicated explicitly, but the constraint of 
fixed end-positions qA and q B is left implicit. The latter convention for the 
end-positions is also followed below in (2.2), (2.3), etc. From the Action Prin- 
ciple (2.1) one can derive Lagrange's equation(s) for the trajectory, which we 
shall not do in this review; see for example [3J El- 



(SS) T = 0, 



(2.1) 
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The reformulated (or General) Maupertuis Least Action Principle (GMP) 
states that, for a true trajectory Maupertuis' Action W is stationary on 
trajectories with fixed end-positions qA and qs and fixed mean energy E: 

{5W) B = 0. (2.2) 
Here Maupertuis' Action W is given by 

IB t B 

W = J pdq = J 2Kdt, 

QA t A 

where p = dL / dq is the canonical momentum, and in general pdq stands for 
P\dq\ + ... + Pfdqf . The second form for W is valid for normal systems, 
for which L = K — V and K is quadratic in the q's, where K is the kinetic 
energy and V the potential energy. The mean energy E is the time average 
of the Hamiltonian H(q,p) over the trial trajectory of duration T, 

1 tB 

E = f J dtH(q(t),p(t)), 
t A 

from the initial position q^ to the final position qs which are common to all 
trial trajectories. Note that E is fixed in (2.2) but T is not, the reverse of 
the situation in (2.1). We use the notation E rather than H because in the 
early sections we restrict [TSj ourselves to dynamical variables where the 
Hamiltonian is equal to the energy K + V; the general case is discussed in 
note [12] and later sections. 

The Reciprocal Maupertuis Principle (RMP), or principle of least (sta- 
tionary) mean energy, is 

(5E) W = 0. (2.3) 

For stationary (i.e. steady-state or bound) motions (periodic, quasiperiodic, 
chaotic) (2.3) is the classical limit of the Schrodinger Variational Princi- 
ple of Quantum Mechanics: (5 < H >) n = 0, (see also sec. 7). The connection 
is clear intuitively: for bound motions, in the large quantum number limit 
the expectation value of the Hamiltonian becomes the classical mean energy 
and the quantum number n becomes proportional to the classical action W 
over the motion. This connection enables the RMP to lend itself naturally 
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to semiclassical applications (sec. 2.3). For periodic and quasiperiodic mo- 
tions the Reciprocal Maupertuis Principle (2.3) is equivalent [12] to Percival's 
principle for invariant tori [TH |2*U] . The RMP is, however, more general and 
is valid for chaotic motions, scattering orbits, arbitrary segments of trajecto- 
ries, etc. Percival's principle has been derived by Klein and co-workers from 
matrix mechanics [21] in the classical limit. Hamilton's Principle also has a 
reciprocal version ^2] - see next section. For a summary of reciprocity in 
variational calculus, see Appendix I. 

The textbook Euler-Lagrange version of Maupertuis' Principle differs 
from (2.2); the constraint of fixed mean energy E is replaced by one of fixed 
energy E, 

(6W) E = 0. (2.4) 

This is not erroneous, since a true trajectory does indeed conserve energy, but 
it is the source of the inconveniences referred to by Jacobi ^T]. Because, as 
we show in the next section, the GMP is equivalent to the HP, and therefore 
to the equations of motion from which energy conservation follows, it will 
become clear that energy conservation is a consequence of the GMP (2.2), 
rather than an assumption as in the original MP (2.4). 

It is well known that conservation of energy is a consequence of symmetry 
under time translation (Noether's theorem), either via the Lagrangian and 
equations of motion [221 > or from the action and Hamilton's Principle [2*3*] . 
Similarly energy conservation can be derived directly from the GMP [T2"] . 

In the early sections we assume L = K—V and that V is time- independent 
(conservative). Later these restrictions are relaxed (see sections 5 and 6). 

2.2 Derivation of Action Principles 

Within Mechanics one can show that different Action Principles such as (2.1), 
(2.2) and (2.3) are equivalent to each other. This we sketch below (see ^2] 
for more details). Alternatively, one can show jUEHEl from Hamilton's Prin- 
ciple (2.1) that the Lagrange equations of motion, in some coordinates are 
essentially Newton's equations in the same set of coordinates. Reference [5] 
shows that (2.4) is equivalent to Newton's equations, quoting an argument 
due to Lagrange, and in ^2] we derive (2.3) from similar arguments. On a 
more lofty level one can derive the Action Principles from Quantum Mechan- 
ics in the classical limit. Dirac and Feynman [21] have derived Hamilton's 
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Principle from a path integral formulation of Quantum Mechanics, and sim- 
ilarly the RMP (2.3) can be derived from wave mechanics (see ^2])- Further 
connections between quantum and classical variational principles are given 
in sec. 7. 

The equivalence of Hamilton's Principle (2.1) with the General Mau- 
pertuis Principle (2.2) is worth spelling out in a little more detail. (The 
equivalence of the MP with the GMP is discussed in Appendix II). The 
starting point of the derivation is the Legendre transformation between the 
Lagrangian of the system and the corresponding Hamiltonian, 

H(q,p) = J2Pn<in- L(q,q). (2.5) 

n 

Here n = 1,2,...,/ runs over the various degrees of freedom, and q stands 
for (qi, qf) and p for {pi, ...,Pf). Upon integration over time along a trial 
trajectory starting at A = (g^, 0) and ending at B = (q B , T) (where we could 
have qs = qA for a closed trajectory), and with the definitions given above 
we obtain 

ET — W — S, (2.6) 

where T is the duration along the trial trajectory. Thus W and S are also 
related by a Legendre transformation. If we now vary the trial trajectory 
and the duration T, keeping and qs fixed, and compute the variation of 
the different terms in (2.6), we obtain 

6S + E6T= 5W-T8E, (2.7) 

which is the relation between the different variations. For conservative sys- 
tems, near a true trajectory Hamilton's Principle is equivalent to the vanish- 
ing of the left hand side. In fact a vanishing left hand side is the unconstrained 
form j2H] of Hamilton's Principle (the UHP - see fo r proof): 

5S + E5T = 0, (2.8) 

where we have used the fact E = E on a true trajectory for conservative 
systems. Eq (2.8) has the unconstrained form S(S + AT) = 0, with A = E 
(the energy of the true trajectory) a constant Lagrange multiplier. For fixed 
T (i.e. 5T = ), we recover the HP (8S)t = 0, and for fixed S (i.e. SS = 0) 
we obtain the Reciprocal Hamilton Principle (RHP) 
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(ST) S = 0. 



(2.9) 



The RHP (principle of least time) is discussed in detail in ref.|12j. 

The right hand side of (2.7) (which therefore must also vanish near a true 
trajectory) gives the unconstrained form of Maupertuis' Principle (UMP) for 
conservative systems: 

SW-T6E = 0. (2.10) 

This has the unconstrained form 5(W + XE) = 0, with A = —T ( the duration 
of the true trajectory) a constant Lagrange multiplier. For fixed E (2.10) 
gives the GMP of equation (2.2) while for fixed W it gives the reciprocal, the 
RMP of equation (2.3). An alternative derivation of the UMP (2.10) is given 
in ref. [12] with traditional variational procedures, employing the general first 
variation theorem of variational calculus, which includes end-point variations 

[131211. 

Equation (2.7) connecting the different variations is reminiscent of equa- 
tions of thermodynamics, where two conditions are derived from a single 
equation. There are many other analogies with thermodynamics, e.g. two 
actions in mechanics analogous to two free energies in thermodynamics, with 
a Legendre transform relation in each case [2E], adiabatic processes and re- 
versible processes in both mechanics and thermodynamics, and reciprocal 
variational principles in both mechanics and thermodynamics (see Appendix 
I). Beginning with Helmholtz jHj, many of these analogies have been explored 
by various authors |2*7j . 

As we have stressed, we hold the end-positions qA and qs fixed in all the 
variational principles discussed above. It is possible to relax these constraints 
by generalizing further the UHP and UMP (see notes [22], [26J, and ;124J), 
but we shall not need these generalizations for the applications we discuss. 

The derivation we have sketched shows that the four principles (2.1), 
(2.2), (2.3), (2.9) are equivalent to each other. Mathematically, the Hamil- 
ton Principles can be regarded as Legendre transformations of the Mauper- 
tuis Principles, since the Legendre transformation relation (2.6) allows us to 
change independent variables from E to T. Thus in Classical Mechanics we 
have a set of four variational principles that are symmetric under Legendre 
and reciprocal transformations: 

(6E) W = rec ^ ty (6W) M = Le ^ re (5S) T = rec ^f y (5T) S = 0. 
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After a short interlude on action principles in semiclassical mechanics we 
proceed to discuss some simple examples. 

2.3 Semiclassical Mechanics from Action Principles 

The RMP ( 2.3) is the basis of a very simple semiclassical quantization 
method. As we shall see, the solution of (2.3) yields, for the true trajec- 
tory or approximations to it, £ as a function of the action W, E(W). For 
bound motions we then assume the standard Einstein-Brillouin-Keller (EBK) 
[2*%] quantization rule for the action over a cycle 

W = {n + a)h, (2.11) 

where n = 0,1, 2, ... , a is essentially the Morse-Maslov index ( e.g. a = 1/2 
for a harmonic oscillator), and h = 2nh is Planck's constant. (EBK (or torus) 
quantization is the generalization of Bohr-Sommerfeld- Wilson quantization 
from separable to arbitrary integrable systems). Thus we have expressed E 
as a function of the quantum number n, E — E n . For multidimensional 
systems, (2.11) is applied to each action Wi in E(W±, Wi, ...), so that we 
get E n i )n2j _, as a function of the quantum numbers rij. This method has 
been found to be reasonably accurate even for nonintegrable systems, where 
strictly speaking the good actions Wi do not exist. The UMP (2.10) can also 
be used to obtain semiclassical expressions for E n i jT a,..- Examples are given 
later in sections 4.1 - 4.3. 

3 Practical Use of Variational Principles. Ped- 
agogical Examples 

Here we discuss how to use variational principles to solve practical problems. 
We use the so-called direct method j2Hj of the calculus of variations (e.g. 
Rayleigh-Ritz), which operates directly with the variational principle and 
makes no use of the associated Euler-Lagrange differential equation. His- 
torically [HO], the method originates with Euler, Hamilton, Rayleigh, Ritz 
and others (see [121 12H] for some early references). We assume a trial solu- 
tion, which contains one or more adjustable parameters a,. We then use the 
variational principle to optimize the choice of the dj. Generally speaking, 
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the more parameters Oj, the better the solution. The direct method is well 
known in quantum mechanics, but is also very useful in classical mechanics, 
as we shall see. We give a few simple examples in this section, mainly for 
bound motions; further simple examples, including for scattering problems, 
are given in refs. [T2*l ITB*] and in the references therein . 

The simplest example is the free motion of a particle of unit mass on 
the line segment x — to x — 1. With the traditional Maupertuis principle 
(2.4) once we fix the energy E = v 2 /2 there is essentially nothing left to 
vary. With the GMP (2.2) or its reciprocal (2.3) we can choose a variety of 
trial trajectories. For example we can take velocity v\ on the sub-segment 
< x < 1/2 and velocity t>2 on the sub-segment 1/2 < x < 1. The action 
W is (t>i + 1>2)/2 and the mean energy E is v 1^2/2, where we can vary V\ to 
change W while keeping the mean energy fixed. We then find that V\ — v%, a 
simple example of energy conservation for this particular trial trajectory. By 
iteration we can now see that energy is conserved for the whole trajectory 
on the segment (0, 1), as a consequence of the General Maupertuis Principle 
(2.2). 

3.1 The Quartic Oscillator 

An instructive example is the one-dimensional (ID) quartic oscillator, with 
Hamiltonian 

v 2 1 

# = — + -Cx\ (3.1) 
2m 4 y ' 

Using a harmonic oscillator trial trajectory with 

x{t) = A sin ut, p{t) = mAu cos ut (3.2) 
we find for a complete cycle with period T = 2n/u) 

4-7T S2-K z m z uj z 
W = TrujmA 2 . 

From (3.3), with the use of the RMP (2.3) we can obtain the "best" 
frequency u> solving (dE/duj)w = 0, which can be substituted into (3.3) to 
obtain E 

r 3CW U/3 771 1 



uJ 



Airm 2 



E = =-{C/m 2 ) 1 l\?,W/A'Kfl z (3.4) 
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and the period T function of mean energy: 

T(E) = 27r[m 2 /(2C£)] 1 / 4 . (3.5) 

This variational estimate for the period is in error by 0.75% when compared 
with the exact result, as noted in reference ^3j- Systematic improvements 
can be obtained by including terms _Bsin3a;t, Csin5ut, etc., in the trial 
trajectory x(t); see [T3| for a related example. 

We have used the RMP (6E)w = here since (3.3) made it very conve- 
nient to do so. The GMP (6W)e — and the HP are also viable for this 
problem, but the original MP (8W)e — is not, since the constraint of fixed 
E leaves essentially no freedom for variation for ID problems |31| . One can 
get around this by relaxing the constraint via a Lagrange multiplier, but this 
is then equivalent to the GMP (see App. II) or the UMP. In 2D etc., the con- 
straint of fixed E does allow some freedom for variation, but the constraint 
is very cumbersome and the GMP is always much more convenient. 

The use of the action principle (2.3) to obtain the mean energy as a 
function of the action (3.4) or the period as a function of mean energy (3.5), 
with harmonic oscillator trial trajectories, is very similar to the use of the 
variational principle in quantum mechanics to estimate the energies of states 
of the quartic oscillator with harmonic oscillator trial wave-functions. We 
show this explicitly for the eigenstates of the Hamiltonian (3.1). In both 
cases an optimum frequency for the trial harmonic case is sought and the 
final errors are also similar. We start by computing the quantum expectation 
value of (3.1) using a harmonic oscillator state with trial frequency uj, 

< n\H\n >= \{n + \)%UJ + ^( — ) V + n + ^), (3.6) 
2 2 8 mw I 

and find the optimum frequency uj by demanding the vanishing of the deriva- 
tive of < H > with respect to the frequency. We substitute this frequency 
in (3.6) to estimate the energy of excited states 

En = \ W 4/3 (^) 1/3 (- 2 + n + \yl\n + \f'\ (3.7) 

Taking the asymptotic limit of (3.7) for n ^> 1 and replacing 2nhn by W 
we obtain the second equation (3.4). This is not as trivial as it seems, since 
(3.7) was obtained from the quantum variational principle (for an arbitrary 
state [TE]). while (3.4) was obtained from a classical action principle. 



13 



This example illustrates neatly that the two variational principles are 
related: the Reciprocal Maupertuis Principle is the classical limit of the 
Schrodinger Variational Principle (see sec. 7 for further discussion). The ex- 
ample illustrates also that the classical action principle can be used as a 
method of approximation, and semiclassically to estimate E n as described in 
sec. 2. 3. According to (2.11), we quantize E(W) in (3.4) semiclassically by 
the replacement W — > (n + as just discussed, this leads to agreement 
asymptotically ( n — > oo) with the quantum result (3.7). Reference (TH] 
discusses the similar relation between quantum and classical variations for 
a different one- dimensional system, with a linear potential. The application 
of classical action principles to other simple systems is discussed in reference 
13j. The application of quantum variations for excited states has a long 
history: the earliest reference is McWeeny and Coulson, who applied it to 
the quartic oscillator described above [32] more than fifty years ago (for later 
references see [To]). 



3.2 The Spherical Pendulum. Precession of Elliptical 
Orbits 

For a pendulum with two degrees of freedom (9, (ft), there is an additional 
term in the kinetic energy compared to that of a plane pendulum, i.e. 

K = ^mL 2 {6 2 + sin 2 6<p 2 ) , V = mgL(l — cos 8) , (3.8) 

where 9 is the polar angle (measured from the downward vertical axis) and 
<p is the azimuthal angle, m is the bob mass, L the length, and g the gravi- 
tational acceleration. 

Using the axial symmetry around the z (vertical) axis, we introduce the 
cylindrical coordinate p = (x 2 + y 2 ) 1 ^ 2 = Lsin9, and the angular momentum 
component l z = mp 2 (f), which is a constant of the motion. We also expand 
cos#, keeping up to quartic terms, which are 0(L~ 2 ). This gives 

K = i v(i + £) + JL , v = i™ y (i + \£) , (3.9) 

where u 2 = g/L. With two degrees of freedom and two constants of the 
motion (l z and the energy), the system is integrable. 

If the quartic terms are neglected in (3.9) by assuming L — >• oo, the (x, y) 
or (p, 0) motion is that of a 2D isotropic harmonic oscillator with frequency 
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uj . In particular, an elliptical orbit, with semiaxes a and 6, is fixed in 
space, and the motion is periodic. If the quartic terms are not negligible, 
the oscillator frequency shifts from ujq to a lower value ui, as in the plane 
pendulum case, where 

( a 2 + b 2 \ . , 

and the orbit precesses in the prograde sense with frequency Q, where 

n = ~^. (3-11) 

The motion is now in general quasiperiodic. The precession is one of the 
major sources of error to be eliminated in the design of Foucault pendula [33] . 
(The Foucault pendulum is discussed in Sec. 5.2.2.) Because the precession 
rate is 0(L~ 2 ), this precession is particularly relevant for recently developed 
short Foucault pendula, which are less than 1 m in length. Equations (3.10) 
and (3.11) have been derived by perturbation theory [34J. 

To derive (3.10) and (3.11) variationally ^3], we first choose a trial tra- 
jectory; we take an ellipse which is closed in the coordinate system x',y' 
which rotates with respect to the fixed x, y system with the (unknown) pre- 
cession frequency Q. In the x',y' system we therefore have x' = acosut and 
y' = b sin ut, where u may differ from uq. Because the x', y' and x, y systems 
differ by a rotation, we have p 2 = x 2 + y 2 = x' + y' , and therefore our trial 
trajectory in the fixed coordinate system is 

p 2 = a 2 cos 2 tut + b 2 sin 2 cut . (3. 12) 

To apply the Maupertuis principles, we calculate the action W and mean 
energy E over one period T = 27r/u of the rotating ellipse, 

W = T{2K), E=(K + V), (3.13) 

where (...) denotes a time average over one period, and K and V are given by 
(3.9). The averages of the terms in (3.9) are easily calculated using (3.12), 
which gives 

(P 2 )= 1 -(a 2 + b 2 ), (p 2 ) = y(a-6) 2 , <i> = i, 
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(3.14) 

<p 2 p 2 > = ^(a 2 -& 2 ) 2 , <p 4 ) = ^(a 2 -& 2 ) 2 + a 2 & 2 . 



The precession frequency f2 can be related to the angular momentum 
using oj + VL = ((f)) = (l z /mp 2 ), and (1/p 2 ) = l/ab from (3.14), giving 

l z = mab(uj + £1) = mabuj ^1 H ^ . (3.15) 




Figure 1: The precessing orbit. We choose positions 1 and 2 as the fixed end 
positions to be used in the action principle, (figure adapted from ref.Q) 

As initial and final positions 1 and 2 (see Fig.l), we take two successive 
aphelia at x' = a, i.e. (pi,<j>i) = (a, 0) and (p2,4>2) = (a, 4>i)i where 02 = 
(u + Q)T = 27r(l + f2/co>).The positions 1 and 2 and hence a and VL/uj are 
thus fixed, which leaves b and to available as variational parameters. Note 
that (3.15) then implies that l z is varied with the trial trajectory (although 
it is a constant of the motion on the individual trial trajectories). We find it 
easiest to use the unconstrained version of the Maupertuis principles [i.e. the 
UMP (2.10)], 5E = T~ l 5W, where T = 2-k/uj is the time to reach position 2 
from position 1. Thus we use 
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where, choosing units such that m = 1, L = 1, coq = 1, we have from (3.13) 
to (3.15) 

E=^(a-6) 2 + ^(a 2 -6 2 ) 2 + y(l + §) 2 ^ 

+ 1 -(a> + b*) + ^-br+ 1 -a%\ (3.17) 

and 

— = -(a - 6) 2 + -(a 2 - 6 2 ) 2 + -(1 + -fab . (3.18) 



The first of the variational equations (3.16) is satisfied identically by 
(3.17) and (3.18), and gives no useful information. From the second, setting 
uj = ujq + A = 1 + A, where A is the frequency shift, and retaining terms only 
to second order in a and b, i.e., 0(L~ 2 ) in general units, we get the condition 

When a = 0, (3.19) gives A = — 6 2 /16 as it should [recall the plane 
pendulum result ^21 EI]]- By symmetry, A must reduce to — a 2 /16 when 
6 = 0. This requires that Q = (3/8)a6 in (3.19), and from this result we 
then get A = —(a 2 + 6 2 )/16 from (3.19). When we restore the dimensional 
quantities ujq and L, these results agree with (3.10) and (3.11), so that with 
a harmonic oscillator trial trajectory, the variational method gives results 
correct to 0(L~ 2 ). 

Hamilton's Principle can also be applied to this problem. Here, we need 
the Hamilton action S = f$ Ldt = T(K — V), where T (and hence uj) is 
now fixed, in addition to the positions 1 and 2 (and hence a and Vt/uj). This 
leaves b available as a variational parameter, and setting dS/ db = leads to 
the same results as above. 



4 More Complicated Examples. Research Prob- 
lems 

More complicated examples are easy to find, and are especially interesting 
in multidimensional systems which exhibit chaotic motion. Such systems 
are nonintegrable, and have fewer constants of the motion than the number 
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of degrees of freedom. We illustrate the use of classical and semiclassical 
variational principles together with the application of quantum variational 
principles to excited states. 

4.1 The x 2 y 2 oscillator 

The 2D x 2 y 2 oscillator is a simple example of this class, which has been 
studied in the literature for some time (see, e.g. jSSl EH])- There is only 
one constant of the motion, the energy. Classically, most trajectories in this 
system are chaotic [37]. Fig.2 shows the contours of the x 2 y 2 potential, and 
Fig. 3 shows a typical (chaotic) classical trajectory. 

6.0 

3.0 
* 0.0 
-3:0 
-6.0 

x 




Figure 2: Contours of the potential V(x, y) = x 2 y 2 /2. The darkened contour 
corresponds to V(x, y) = 1. (from ref.[3S]) 

Refs. [I2| HSj give classical and semiclassical discussions of the motion us- 
ing the RMP. Here we discuss the quantum states using variational methods 
for the excited states, as done above for a simpler example. Comparison with 
the semiclassical results is given later. The Hamiltonian is 

H = K x + K y + V(x, y) = P ^±i + ^f. (4.1) 
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z 



Figure 3: A typical trajectory of the Hamiltonian (4.1). The motion starts 
near the origin, then proceeds into the right-hand arm of the potential along 
the positive x-axis. The trajectory eventually reverses its x-direction and re- 
turns to the vicinity of the origin before entering another arm of the potential, 
(from ref . |35j ) 

We choose units where m = 1 and C = 1 , where Cx 2 y 2 /2 is the potential. 
Below we also choose h — 1. We take as an approximation a separable trial 
wave-function ip nx ,n v (x,y) 

^n x ,n y (x,y) =ip nx (uj x ;x)ijj ny (uj y ;y), (4.2) 

where the wave-function ip nx (uj x ; x) is an eigenstate of a harmonic oscillator 
with frequency uj x : 

p 2 -\- uP'x 2 1 
x 2 x ^n x {u x ;x) = {n x + -)u x ip nx {u x -x). (4.3) 

The frequencies u x and u y are variational parameters, which will differ 
from state to state, and as a result the set (4.2) is nonorthogonal. The com- 
putation of matrix elements and expectation values with the wave functions 
(4.3) is rather easy, 
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< >= -(n + < n\x 2 \n >= — , (4.4) 

2 2 2ui 

and therefore the expectation value of the Hamiltonian (4.1) with the wave- 
functions (4.2) is 

11 1 (n x + -)(n + -) 
< H >=< n x ,n y \H\n x ,n y >= -[(n x + -)Ld x + (n y + -)ujy^ 1 

2 2 2 UJ x UJy 

(4-5) 

We then find the best frequencies oo x and oo y by setting to zero the partial 
derivatives of < H > with respect to the two u;'s to obtain the frequency 
conditions 



1 ^ 2 / 1 

-) , UxUy = {n x + - 

from which we obtain 



u 2 x Uy = (n y + -) , (j x Uy = (n x + -), (4.6) 



(n x + ^)uj x = (ny + ^)uj y . (4.7) 

Eqs. (4.7) and (4.3) show that at the variational minimum the ratio of the 
two frequencies is rational and the mean energy of the re-oscillator is the 
same as the mean energy of the y-oscillator. It follows immediately that the 
mean kinetic energies of the two oscillators are also equal to each other at 
the variational minimum: < K x >=< K y >. From (4.6) we also find 



1, , 1 



UO x U). 



x^y 



2 

which gives the separate frequencies 



[K + o)K + o)] 1/3 > (4-8) 



<"* = K + \)' l/ \n y + l -f\ u y = (n x + l -fl\n y + (4.9) 

and the variational estimate for the energy: 

3 11 

E nx) n y =< n x , n y \H\n x , n y >= -[{n x + -){n y + -)] 2/3 . (4.10) 

Therefore the variational states are labelled by the two integers n x and n y 
and the energy at high quantum numbers increases as the two thirds power 
of their product. 
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The estimate (4.10) can also be obtained semiclassically by using the 
RMP to find E(W X , W y ), and then quantizing the actions Wi with the EBK 
relation (2.11); see ref.^3] for a detailed discussion. It is to be noted that 
this method yields the quite reasonable first approximation (4.10) for E nxiUy 
despite the fact that strictly speaking the good actions W x and W y do not 
exist for the x 2 y 2 oscillator. In essence, one is approximating a chaotic tra- 
jectory by a quasiperiodic one, and the RMP picks out the best such trial 
trajectory |3~Hj . 

The estimate (4.10) differs only in the leading coefficient , 3/2 vs. 1.405, 
from a semiclassical adiabatic formula (SCA) obtained by Martens et al [33] . 
The difference between the two estimates is about 6.5%. As noted by Martens 
et al [35] formula (4.10) gives rise to degeneracies. For example (n x ,n y ) 
= (7,0), (0,7), (2,1), (1,2) give the same energies - a four-fold degeneracy. 
However, numerical estimates of the first 50 energy levels, obtained in |3*3] . 
do not show such high degeneracies. Also, because the symmetry group of 
the Hamiltonian C& v (a fourfold axis of rotation, and four reflection planes) 
has only one- and two-dimensional irreducible representations j3H], we expect 
at most two-fold degeneracy. We shall see that the basis (4.2) permits an 
evaluation of the splitting of these spurious degeneracies. The pattern of 
degeneracies is easy to understand if the bracket is rewritten 

[(2n x + l)(2n y + l)/4]=N/4, 

where the number N (the "principal" quantum number), is an odd integer. 
The degeneracy pattern is determined by the decomposition of N into prime 
factors. For example (n x , n y ) = (2,l) corresponds to N =15 which can be 
factored as 15 • 1, 1 ■ 15, 3-5, 5-3 and these products correspond to (7,0) 
etc. If N is a prime number there are only two states, corresponding to N ■ 1 
and 1 ■ N, while every nontrivial factorization of a nonprime N allows further 
states. The integer N belongs to one of two possible classes: N = 4k + 1 
or iV = 4k — 1. For the case 4k — 1 (e.g. iV =15) one of n x ,n y is even and 
the other is odd. This class of states EO and OE (of 35J) corresponds to 
the two dimensional irreducible representation of C± v . For N = 4k + 1 ( 
e.g. N—21) the two quantum numbers (n x ,n y ) are both even or both odd 
(e.g. (n x , n y ) = (10,0) or (0,10) or (1,3) or (3,1) for N—21, which is therefore 
four-fold degenerate). This degeneracy is split into states of symmetry EEE, 
EEO, OOE and OOO, in the notation of Martens et al [35] for the one 
dimensional irreducible representations. In this notation the first two letters 
give respectively the parity ( Even or Odd) of n x and n y (or equivalently the 
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behaviour of the wave-function under the reflections (x — > —x, y — > — y), 
while the third letter stands for the behavior under the interchange of x with 
We discuss later the splitting of degeneracies resulting from the expression 
(4.10). At this point we wish only to note that the degeneracy grows very 
slowly with N, like lnlniV, and even for N quite large where quasiclassical 
behaviour is expected (say iV ~ 105 ) the typical number of degenerate states 
is quite modest, 4-6. 

The main points of discussion arising from the result (4.10) are the split- 
ting of the degenerate levels, perturbative corrections to the variational en- 
ergies , and the density of states implied by this formula. For the first fifty 
low lying states one has the accurate numerical estimates of Martens et al 
35j with which to compare. 



4.1.1 Perturbative Shift of Variational Results 

We can improve the variational result (4.10) if we rewrite the Hamiltonian 
H (4.1) in the form 

u 2 x 2 + u 2 y 2 x 2 y 2 - Jix 2 - u 2 y 2 
H = H + V = K + — yJ- + ^. x - V—. (4.11) 

Here, the perturbation V is the last term above, and its effect is already 
included in first order, in the variational estimate (4.10). Therefore the first 
correction comes in second order in V. We must also remember that that 
the frequencies (u x ,u y ) are appropriate to a fixed pair of quantum numbers 
(n x ,n y ) , and change if we go to a different pair. The simplest way to do 
perturbation theory in these circumstances is to estimate corrections to each 
state (n x ,n y ) by considering a set of orthogonal intermediate states (n x ,n y ) 
all with the same frequency (u x ,w y ) appropriate to the state (n x ,n y ) which 
we consider. Therefore we take the second order correction 



< n' x ,n y \V\n x ,n y > 
[« - n x )u x + (n - n y )u y ] 



A4 2 > = - v S' n ry\ (4.12) 



Due to the form of V (quadratic in the coordinates x, y), we must have 
n x = n x or n x = n x + 2 or n x = n x — 2 for nonzero contributions and 
similarly for n y . It is an interesting property of the variational solution that 
only intermediate states with both n' x different from n x and n' y different from 
n y give nonvanishing contributions, because of the frequency conditions. This 
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can be seen from an evaluation of matrix elements with only one quantum 
number different, for example n x ^ n x , but n y = n y . Then one finds 

< n' x ,n y \V\n x ,n y >=< n x \ — \n x > [< n y \y 2 \n y > (4-12') 

and the bracket on the right hand side vanishes because of (4.6). Only the 
perturbation x 2 y 2 /2 contributes to nonvanishing matrix elements. There are 
only four intermediate states which contribute: (+, +), (+, — ), (— , +), (— , — ), 
where (+, — ) = (n x + 2,n y — 2), etc. Therefore (4.12) becomes after a little 
computation 

AE£U = -^ Ei nL v , (4-13) 

and 

E$, ny = —{{n x + 1/2)K + 1/2)] 2 / 3 , (4.14) 

where AE^ denotes a second-order shift and E^ denotes the energy up to 
the second-order. 

Therefore, in second-order perturbation theory, the coefficient 3/2 is di- 
minished by about 4% to 1.437, closer to the adiabatic result 1.405 of Martens 
et al [35] . The remaining difference is about 2.5%. The degeneracies of (4.14) 
are still the same as those of (4.10). If we apply perturbation theory to the 
ground state (0, 0) we find 



E 00 = E ( w> + AE { Z > + AE$ + AE$ = 0.5953 - 0.0248 + 0.0496-0. 1465 + ... 

(4.15) 

The sequence (4.15) indicates that the perturbation expansion does not con- 
verge, and it is a good idea to stop at second order. The accurate numerical 
result for £"00 (from [33] ) is 0.5541, which shows that stopping at second order 
gives an error of about 3% for the ground state. 

Some examples of the estimates E^\E^ 2 \ and Esc a are given in Table 
1, and compared with numerical estimates from jSH], for a set of states which 
are also relevant to the next section. 
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Table 1 



N 


[Tlx j Tly) 


E eX act 


E(°) 

var 


E (2) 


E S plit 


Esca 


1 


(0,0) 


0.5541 


0.595 


0.570 


na 


0.5577 


5 


(2,0) 


1.7575 


1.7406 


1.6681 


1.7471 


1.6308 


9 


(4,0) 


2.4920 


2.5756 


2.4683 


2.5377 


2.4132 


13 


(6,0) 


3.1180 


3.2911 


3.1540 


3.2171 


3.0836 


17 


(8,0) 


3.6909 


3.9357 


3.7717 


3.8309 


3.6874 


21 


(10,0) 


4.2504 


4.5310 


4.3422 


4.3973 


4.2453 


25 


(12,0) 


5.0490 


5.0895 


4.8775 


4.9297 


4.7685 


29 


(14,0) 


5.3327 


5.6189 


5.3847 


5.4347 


5.2645 



EEE states at low excitation. The exact values and Esca are taken 
from [33]; the variational estimates and second order results correspond to 
formulae (4.10) and (4.14), while the E sp n t are obtained following section 
4.1.2 

4.1.2 Splitting of Degenerate States 

The basis (4.2) is also quite convenient to obtain the splitting of degenerate 
states. We shall deal with the simplest example with a degeneracy: the pair of 
states (2, 0) and (0, 2). Both these states have the same energy (3/2)(5/4) 2 / 3 
in the variational approximation. The wave-functions of the two states are 
nonorthogonal, because the frequencies in the x and y oscillators differ. We 
define the overlap integral A 

A =< 20|02 >= | J dxMuj 2 ;x)M^;x)\ 2 = ^(\f'\ (4.16) 

where the two wave functions have different frequencies ( lo-i = (1/10) 1 / 3 and 
coo = (25/2) 1//3 ) as denoted. The splitting of two (or more) nonorthogonal 
degenerate states is well known; a simple example is the Heitler-London 
treatment of the hydrogen molecule. One has to take the sum or the difference 
of the two nonorthogonal states , e.g. = N + [ip 02 + 1P20], and similarly for 
The eigenvalues have the form 

E + = (H AA + H AB )/(1 + A), 

and 

E_ = (H AA -H AB )/(l-A), (4.17) 
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where here \A >= |2, > and \B >= |0, 2 >. 

Thus in our case we need the off-diagonal matrix element of the Hamil- 
tonian between the two states (0,2) and (2,0). This can be obtained from 
the overlap integral (4.16) 

2 2 ^) 2 2 1 

< 201X3 + ^102 >= -u 2 A, <2Q\K x + ^-\02>= -cu A. (4.18) 
If we subtract the two equations (4.18) we obtain 

< 20\(u% - c^j)y |02 >= (5lu 2 - u )A. (4.19) 

The right hand side of (4.19) vanishes because of the frequency condition 
(4.7) and therefore the matrix element of x 2 vanishes as well between the 
two variational states. The same is true for matrix elements of y 2 , and x 2 y 2 . 
Therefore only the kinetic energy terms contribute to the off-diagonal matrix 
elements of the Hamiltonian . Using equations (4.18) we find 

< 20IXJ02 >= 12 y° A 9 = -u A. (4.20) 
u$-u4 2 

The same value is obtained for the matrix element of K y , and therefore 
we find 

< 20|# |02 >= u A. (4.21) 
In our case u = (4/3) E 02 , and therefore we have for the two split eigenvalues 

B _ = ^1^^ = 1.8254 (4.22) 

Hence the two degenerate states at E 02 = 1-7406 split into an EEE state at 
1.823 and an EEO state at 1.625. We identify as an EEE state because 
has EEE symmetry, and similarly 1 J r _ has EEO symmetry. When one 
compares with the numerical results of Martens et al [3*3] . the estimates 
(4.22) are qualitatively correct; the error of E + is about +3.7% while the 
error of E_ is +6.4%. One can further apply the factor of 23/24, obtained 
in the previous section, to end with the estimates 1.7471 and 1.556. The 
result E = 1.7471 for the (2,0) EEE state is given in Table I as E sp u t . Other 
low-lying EEE state energies are also given. 

We end this brief excursion into the x 2 y 2 model with the conclusion 
that the quantum variational method for excited states is particularly 
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well adapted for energy estimates in this system. We have used variation- 
perturbation theory here to demonstrate this method. The same results can 
also be obtained purely variationally, by using more elaborate trial states 
than (4.2); e.g., for an EEE state we would try (for n x ,n y even) 

and similarly we would try ip^ n (x,y) for the EEO states. We return to 
purely semiclassical variational methods in Sec. 4.3. For completeness, we 
note that classical variation-perturbation methods have also been devised 

EH- 

4.2 The "Linear Baryon" 

Another system very amenable to classical and quantum variational estimates 
is a system of three equal mass particles constrained to move along a straight 
line |41j . There are constant attractive forces between the particles, or a 
linear potential ( as for ID gravity ), and the particles can move through 
each other. In suitable units the Hamiltonian can be expressed as 

H = ^[Pi+pI+pI] + ^[\ x i - x z\ + \%2 - 2:3 1 + 1^3 - (4.23) 

where Xi are the positions and pi the momenta. If we keep the centre of mass 
at the origin we have the constraints x\ + xi + X3 = and pi + P2 + P3 = 0. 
We can then take internal Jacobi coordinates p, A defined by X\ — £2 = pV% 
and Xi + X2 — 2x3 = Av^6, and the Hamiltonian becomes 

H = \\p 2 p + p\\ + ^[|AV3 - p\ + \\VS + p\ + 2\ P \\. (4.24) 

Taking now polar coordinates R, 9 in the p, X plane, where A = Rcos9 and 
p = RsinO, we obtain 

H = l\p 2 R +Pe/R 2 ] + ^=R(\ sin^| + | sin(^ + 27r/3)| + | sin(^-27r/3)|). (4.25) 
2 y 2 

The equipotentials of this Hamiltonian form regular hexagons. The six 
sides correspond to the six different permutations of the three particles in the 
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baryon; for example, one side corresponds to x\ < %2 < X3, etc. The Hamil- 
tonian (4.25) describes the motion of a particle inside this regular hexagon. 
This system has been studied in the literature [12] under the restriction to 
a single "wedge" (say < 9 < 7r/3) under the name of " wedge billiards 
in a gravitational field" . Recently also the relativistic generalization of the 
Hamiltonian (4.23) has been discussed [43J. 

The classical motion with the Hamiltonian (4.25) can be solved piecewise 
exactly. The trajectory consists of segments of parabolas in each wedge of 
the hexagon. One can visualize the system as a hexagonal "pit" or "funnel" 
in a gravitational field , with a point particle sliding inside the funnel. Since 
the system is not harmonic there are no normal modes, but an infinity of 
trajectories at each energy. Many of these can be simply described in terms 
of two prototypes having high or low "angular momentum" respectively. We 
use quotation marks since angular momentum is not conserved in a hexagonal 
potential. However the average angular momentum of specific orbits is still 
meaningful. The low "angular momentum" trajectory consists in a motion 
along a groove of the hexagonal funnel corresponding to 9 = 0, to the origin 
and continuing along 9 = n, and back. The high "angular momentum" 
trajectory is composed of segments of parabolas joined together to look like 
a circle. Each of these trajectories can sustain transverse oscillations to give 
rise to quasiperiodic motions. Examples are shown in references |41[ 14*3] . 
Eventually for large transverse oscillations the two types of trajectories blurr 
and and chaotic motion ensues. 

With only one constant of the motion (the energy) and two degrees of free- 
dom, the system is nonintegrable, and possesses chaotic trajectories. These 
and the nonchaotic ones are piecewise parabolic, but it may be advantageous 
to represent them by simpler trial functions. For illustration we describe a 
low angular momentum trajectory 9 = pg = when (4.25) becomes 

H = f + (4.26, 

Consider the trajectory with E — 1. The problem can be solved exactly to 
give a period (for the whole cycle) of T exact = 4.6188, while a variational 
estimate with the RMP using a harmonic oscillator trial trajectory (with the 
same mean energy E = 1 as the exact trajectory) gives T variat i ona i = 4.6526, 
an overestimate of 0.73%. 

Quantum variational methods are quite useful here too since we can- 
not solve the quantum problem exactly, either for the ground state or ex- 
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Figure 4: The quantum probability distributions for four excited states of 
the hexagon pit. The notation nl2m0 means the state n = 12, m = 0. The 
color scale of relative probabilities is given at the bottom. The graphs have 
been rescaled to make the hexagon, the classical boundaries of the orbits for 
the given E, the same size, (from ref. [H]) 

cited states. To describe the symmetry species of the symmetry group Cq v 
of the quantum Hamiltonian we can use eigenf unctions of the 2D isotropic 
harmonic oscillator 03]. In polar coordinates these are products of an ex- 
ponential exp(— uR 2 /2) times a function of 9 ( sin(m#) or cos(m#)) times 
(cuR 2 ) m/2 times a La guerre polynomial in uR 2 times a normalization factor. 
Although the azimutal quantum number m is not conserved with the Hamil- 
tonian (4.25) it is conserved modulo 6 due to the hexagonal symmetry of the 
problem. We can as a first approximation minimize the diagonal elements 
of the Hamiltonian (4.26) in this basis with respect to the frequency and 
obtain first order estimates for the energies of states which are labelled by a 
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principal quantum number n and m(modulo 6). For example the energy of 
the ground state (n,m) = (0,0) is overestimated by about 0.39% in this basis 
while energy of the state (n,m) = (15,0) is overestimated by about 3%. "Ex- 
act" energies can be obtained by diagonalizing the Hamiltonian in a large 
basis of harmonic oscillator states and minimizing eigenvalues with respect 
to uj. Another variational procedure uses products of Airy functions along 
the groove and perpendicular to it, and takes the scale factors of the two 
Airy functions as variational parameters. For the excited state (15,0) the er- 
ror is 0.014% (the "exact" energy is 16.34869 while the Airy approximation 
gives 16.3511). More examples can be found in Table 1 of reference jJJJ . The 
probability densities of a few states are shown in Fig. 4 where one can see one 
low angular momentum state (12,0) and one high angular momentum state 
(2,18) and two states which presumably correspond to chaotic trajectories. 

4.3 He Atom Inside C70 Cage 

Endohedral fullerene complexes such as He @ C 70 have attracted much study 
in recent years, both experimental and theoretical (for reviews see |45| 021 
E7J HH1 EH]). Here X @ C n denotes species X (atom, ion or molecule) encap- 
sulated in a C n cage [5Tj|. Early theoretical work focussed on the prediction 
of infrared and Raman [SHI 134] and neutron scattering [49J spectra due to 
the motion of X = He, Li + , Na + , K + and CO inside the near-spherical Cgo 
cage. The predictions of (33] were subsequently confirmed experimentally 
|5l)l 157)] for the infrared bands of Li + @C6o- Confirmation of the predicted 
Raman bands jSH] of Li + @C6o is tentative [SB]. Later the theoretical studies 
were extended to study the dynamics and spectra of atoms, molecules and 
ions trapped in nonspherical cages, e.g. He @ C70 [HZI EE1 EH] , Ne @ C 70 jEU], 
He @ C 250 (nanotube) [61], Na+ @ C 120 (nanotube) [62], and even (Cgo) n @ 
Coo (nanotube), so-called nanopeapods jHH]- Figs 5 and 6 show Ceo, C70 and 
carbon nanotube cages. 

The C 60 cage, assumed to be rigid, presents a very nearly spherical poten- 
tial [64J for an entrapped atom or ion X and the dynamics has been studied 
by various classical [SH], semiclassical [48J and quantum 0H] methods. The 
potential can be calculated by a variety of semiemperical and ab initio meth- 
ods jlH HHl EH!! El] ■ The dynamics is usually quite regular [64J for X@Ceo- 
We discuss here instead the case He @ C70, where the potential is highly non- 
spherical and gives rise to nonintegrable (chaotic) dynamics of the He atom 
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Figure 5: (a) The soccer ball-shaped molecule Ceo (icosahedral symmetry 
group I/j); (b) The rugby ball-shaped molecule C 70 (symmetry group D 5rf ). 
(from ref. 

jHZl EE1 EH] ■ We estimate the energy levels semiclassically jHHl ESI usm g the 
UMP as described earlier, and compare with the "exact" results obtained by 
brute-force diagonalization of the Hamiltonian jSTj. Classical simulations of 
the motion have also been carried out [H3 EHl EH] , and provide insight as to 
the degree of chaos in the motion of the He atom. The chaos is surprisingly 
small, given that the potential is strongly anharmonic. The small degree of 
chaos implies an approximate third constant of the motion [68] , in addition to 
the energy and ^-component of angular momentum (assuming the potential 
is axial). 

The potential for He in a rigid C70 cage has been estimated semiem- 
perically [HZ], and, to a good approximation, is axially symmetric around 
long z-axis in Fig 5b . With origin at the centre of the cage and x, y axes 
perpendicular to the symmetry z-axis, the potential has been fit to the form 

V(x, y, z) = l -h{x 2 + y 2 ) + h 3 z 2 + C{x 2 + y 2 )z 2 . (4.27) 

The first two terms describe a 3D anisotropic harmonic oscillator, with fre- 
quencies ojW = cj(°) = (fcx/m) 1 / 2 = 84.166 cm -1 and uf> = (/c 3 /m) 1//2 = 
21.934 cm -1 , where m is the He atom mass. We use the standard spec- 
troscopic unit (cm -1 ) for frequency and energy (see (SHj)- The last term 
C(x 2 + y 2 )z 2 in (4.27) is a quartic perturbation, with C = 1.738 cm -1 when 
x,y,z are in units of (h/muj^) 1 ^ 2 . There are also smaller terms in (4.27) 
varying as (x 2 + y 2 ) 2 , z A , etc, which we shall ignore. 

Because of the axial symmetry we change to cylindrical polar coordinates 
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Figure 6: Carbon nanotube of diameter 0.7 nm. This is a quite narrow 
nanotube with the hemispherical endcaps derived from Cqq. The molecule 
depicted is C 250 ; actual nanotubes are much longer [S3], up to Ci )0 oo,ooo an d 
beyond, (from ref [01]) 

(p, (f>, z) and write the Hamiltonian, with potential (4.27), as 



where p = (x 2 +y 2 ) 1//2 , and l z is the z-component of the angular momentum, 
which is a constant of the motion. The cage has been assumed to be rigid 
and fixed in space. To apply the UMP (2.10), we proceed exactly as for 
the simpler spherical pendulum example of sec. 3. 2. When viewed down the 
z axis, the trajectory is roughly a precessing ellipse, and the z-motion is 
oscillatory. We therefore take clS cL trial trajectory 

x = A x cos u x t, y = A y sin u> x t, z = z = A z cos u z t, (4.29) 

where the x ,y axes precess with some (unknown) frequency Vl with respect 
to x, y, and where cu x may differ from a4°), etc. Calculating E and W with 



the trial trajectory and applying the UMP gives to first order in C (see joH 

for the details) Q = 0, the shifted frequencies 



H = -m(p 2 + z 2 ) + 




- + \m^fp 2 + \m^fz 2 + Cp 2 z 2 , (4.28) 




(0) 



CW Z 
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and also, to first-order in C, E as a function of the partial actions Jj = Wi/2ir 



E(I X , I v , I z ) = (4 + J>0°> + I z uf + ^l^; . (4.31) 

The terms of second-order in C in (4.30) and (4.31) are also derived in 
and it is also found that Q = to second order. The simulations show 
that Q is small but nonzero in general, and to obtain Q ^ from the UMP 
requires an improved trial trajectory over (4.29), e.g. one involving higher 
harmonics of uj x and u y [59J. We shall ignore this small effect here. 

To estimate the energy levels semiclassically from (4.31) is very simple: 
according to the EBK quantization rule (2.11), we replace I x by (n x + \)h, 
etc. This gives 

E Wz = (n a +n,+l)^)+K+l)M 0) + CK + % 2 + (0 wo; + ^ +0(C 2 ) 

(4.32) 

where the 0(C 2 ) and 0(C 3 ) terms are written down in ref. jSH]- We compare 
the semiclassical energy levels (4.32) (including the 0(C 2 ) and 0(C 3 ) terms) 
with the fully quantum results; the lowest 613 levels, up to the energy of 
1000 cmT x are obtained by diagonalizing the Hamiltonian in a harmonic 
oscillator basis in [SB]. The harmonic oscillator chosen has a potential given 
by the harmonic part of (4.28). The semiclassical and fully quantum results 
agree to better than 1%. We illustrate the good agreement in Fig 7, which 
shows the levels up to 500 cm -1 , organized according to m, the quantum 
number corresponding to l z . 

For an isotropic 2D harmonic oscillator, the energy depends only on N xy = 
n x + n y . For N xy = 2, for example, we can have (n x , n y ) = (2, 0), (1, 1), (0, 2), 
which are degenerate. The m-values allowed for N xy = 2 are m = —2, 0, +2. 
For general N xy , we have \m\ = N xy , N xy — 2 , N xy — 4 , . . . , 1 or . For our oscil- 
lator (4.28) the different \m\ states should in reality have different energies; 
this is found in the fully quantum calculation, but not in the semiclassical 
one where the energy depends only on N xy = n x + n y . Table 2 illustrates 
more clearly the slight \m\ splittings found quantum mechanically, but not 
semiclassically. 

To obtain the |m| splittings semiclassically requires an improved trial 
trajectory over (4.29), one that gives a nonvanishing precession frequency Q. 
This is discussed in ref jSU]. The simple improvement in the trial trajectory 
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Figure 7: Comparison of energy levels calculated via quantum mechani- 
cal methods (black) and a semiclassical variational technique (grey). (from 
ref.jSHl) 

mentioned earlier (including some higher harmonics) does indeed give a non- 
vanishing Q and hence \m\ splitting. However the trial trajectories employed 
thus far are unable to give accurate splittings, and hence require further 
improvement. 
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Comparison of selected energy levels calculated via quantum and semi- 
classical methods demonstrating the failure of the simple variational calcula- 
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tion to produce level splitting observed in the quantum results, (from ref . [5%] ) 

5 Nonconservative Systems 

In this section we consider nonconservative mechanical systems. For example 
consider a system of particles moving in the time-dependent external field, 
or a system whose parameters depend explicitly on time. Other examples 
include time-dependent constraints, e.g. a bead on a rotating hoop whose an- 
gular velocity is a prescribed function of time, and variable mass systems (for 
which L ^ K—V 'm general |69j). Such systems are not closed and the energy 
is not conserved. The evolution of a nonconservative system is completely 
determined by a time-dependent Lagrangian L(q, q, t) or by a time-dependent 
Hamiltonian H(q,p,t). (We exclude dissipative systems from our discussion. 
We do note that Lagrangians and Hamiltonians have been found for such 
systems in some cases [70J, so that the action principles discussed here would 
apply in these cases. Alternatively [5], dissipative terms can be added to the 
Lagrange equations of motion or to the Hamilton principle.) H and L are 
related to each other by the Legendre transformation (see eq. (2.5)) 

H(q,p,t) = ^p n q n - L(q,q,t). (5.1) 

n 

Locally these lead to the usual Euler-Lagrange or Hamilton equations of 
motion that govern the dynamics of the system. 

5.1 Action Principles 

We consider the global description of nonconservative systems. It is well 
known that the Hamilton Principle 

{SS) T = (5.2) 

is applicable for nonconservative systems as well as for conservative ones. 
Here subscript T reminds us that variations near the actual path are con- 
strained by fixing the travel time T from the initial point A = (q A , 0) to the 
final point B = (q B ,T). Again we have chosen t A = and t B = T, and, as 
before, it is understood here and below that the end-positions q A and q B are 
also fixed. 
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As before, we can relax the constraint of fixed T and rewrite f5.2) |12[ l7T] 
as an unconstrained Hamilton Principle (UHP) 

6S + E(T)6T = 0, (5.3) 

where the Lagrange multiplier E(T) is the true value of the energy at the 
time of arrival at the final point B. (Elsewhere in this paper E(T) is denoted 
variously as E(tB), Eb or Hb] see notes |25J and [IS])- We assume that all 
trial trajectories start at the same time t = and also have the same end 
positions qA and qs- 

Using eq.(2.7), which is derived from (5.1) as before, we transform the 
UHP (5.3) to the form 

T5E -SW- [E(T) - E(T)]5T = 0. (5.4) 

This is the unconstrained Maupertuis Principle (UMP) for nonconservative 
systems, which seems not to have been given in the literature explicitly (it 
is implicit in an Appendix in ^2] - see also |7J3)- It relates three variations 
- the variation of mean energy E, the variation of the action W and the 
variation of the travel time T. The Lagrange multipliers are the true travel 
time T and the difference of energy and mean energy of the true trajectory 
at the final point, [E(T) — E(T)\. The latter can be rewritten as 

T 

E(T) - E(T) = -J dt[t—} =< t— >, (5.5) 
o 

where the bar and brackets < ... > both denote a time average over the 
interval (0,T). The generalized UMP (5.4) can also be derived by traditional 
variational procedures ([121, eq.(A.6)). 

One can derive constrained principles from (5.4) as before, e.g. (5E) WjT = 
0, but these would be complicated to use in practice because of the severe 
constraints on the trial trajectories. The version (SW)^^ = is closest to 
the original Maupertuis Principle. The constrained principle 

Lw + T<t^>}) =0 

is mentioned by Whittaker [73], and is easily derived from (5.4) and (5.5). 
But since (5.4) is also valid without fixing E(0) and E(T) , (5.4) is more 
convenient. 
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Formally the unconstrained Maupertuis principle (5.4) is equivalent to 
the Hamilton principle, but it also appears too complicated to be useful in 
practice. This is not quite true and we demonstrate how the UMP can be 
used for adiabatic perturbations. 



5.2 Adiabatic Invariants and the Hannay Angle 

For simplicity we start with a one-dimensional periodic system that depends 
on a parameter A. Let A vary with time adiabatically, i.e. the change AA <C A 
during a period of oscillation T. It is well known that in this case the energy 
E as well as the period of oscillation T depend slightly on time, but the 
action I over a cycle, where 

1 If 

I = —W(cycle) = — f pdq, (5.6) 

remains constant to an excellent approximation (adiabatic invariant) when 
the parameter A varies. 

Consider what kind of information can be extracted from the UMP (5.4). 
Both the energy E and the period T of an actual motion are functions of the 
action / (recall the action-angle variables): E = E(I), T = T(I). Consider 
two actual trajectories with actions / and I + dl as two trial trajectories 
for the UMP (5.4). Then the variational equation can be rewritten as a 
differential equation 

dE 2tt dE 1 dT , . 

~dl 'T ~ < t ~dt > T~dl (57) 

It is more convenient to rewrite this in terms of the frequency = 
2n/T(I), 

* . dE dE 1 du(I) , . 

^ = !T +<t m > W)v- (5 - 8) 

This equation generalizes the standard relation [S] for integrable conservative 
systems 

dE 

" (/) = III 
to integrable nonconservative systems. 

Up to this point (5.7) and (5.8) are general and valid for any noncon- 
servative system for which a good action can be defined. Consider now the 
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adiabatic regime where the variation of the parameter X(t) over the short pe- 
riod T, AA = A(T) — A(0), is much less than the parameter itself: AA <C A(0). 
In this case 

< t a y» >=< t 8 -^>X >= X < t°*&* >, (5.9, 

where the averaging < tdE(I,X)/dX > is to be performed at a fixed value 
of A ~ A(0) since A varies little over a period T. The action / is conserved 
for any given parameter A, and hence the derivative dE(I,X)/d\ is also a 
constant of the motion. As a result in this approximation we get 

dE(I, A) _ T dE(I,X) _ n dE(I,X) 
< t dX >_ 2 dX ~ u dX ' 

In the same approximation we have for E =< E > in (5.8) 

<E(I,X(t)) >=E(I,\(0)) + 1\^. 

Thus to the first order approximation in A we can rewrite (5.8) as 
, T . dE(I,\(0)) . ...d. 1 dE. dE d . 1 

or finally as 

w(I) = wW + ^. (, 10 ) 

Here w (0) (/,A) = dE(I , X(0))/dI is the zero-order term, and AA = A(T) - 
A(0) ~ X(2n / uj^) is the variation of the parameter A during the given period. 
We get corrections (of the first-order in the nonadiabaticity parameter A) to 
the frequency of oscillation, for adiabatic systems. The result seems very 
simple. It appears that is a continuous function of time t and we 

calculate some average frequency over the period. But actually is not 
always a continuous function of time; it is a global parameter of a (slightly) 
aperiodic system. Thus (5.10) is nontrivial in general. 

To illustrate these results it is instructive to consider a simple exactly 
solvable example. We consider the ID motion of a perfectly elastic ball of 
mass m bouncing between two planes. We suppose that one plane is fixed 
and the other is slowly moving with constant velocity u. We suppose also 
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that the planes are heavy enough so that they are unaffected by collisions 
with the ball. The distance between the two planes L(t) = L + ut is the 
adiabatic parameter. 

Let us define a cycle as the motion of the ball from the fixed plane to 
the moving plane, collision, and the motion of the ball in the backward 
direction up to the collision with the fixed plane. The n-th cycle can be 
described in the following way. At time t — (beginning of the cycle) 
the ball has velocity v n , and the distance between planes is L n . At time 
T n (— >) = L n /(v n — u) the ball collides with the moving plane and after the 
collision it moves to the fixed plane with velocity (v n — 2u) = v n+ i (the 
collision with the fixed plane is elastic and the backward velocity is equal to 
the forward velocity after collision, i.e. f n +i)- The backward motion takes 
time T n (<— ) = (L n v n )/(v n — u)v n+ i. The period of the whole cycle is T n = 
T n {-^) +T n (<— ) = (2L n )/v n+ i. At time t = T n (the initial time for (n+l)-th 
cycle) the distance between the two planes is equal to L n+ i = L n v n /v n+ i. 
This means that L n v n = L n+1 t> n+1 is an exact (adiabatic) invariant. 

It easy to see that the action / is proportional to this invariant, 



2tt 



f P d( l = — L nVn, 



i.e. the action I is an exact (adiabatic) invariant for our special initial con- 
dition [71j . 

In terms of the action / and the adiabatic parameter L the complete set 
of relevant quantities is as follows: 

- 1 ,TcI. ,7rl . 

E n = —(—) 2 -u(—), 

i.e. u) n (I) is linear function of the nonadiabaticity parameter L = u and 
the approximate relation (5.10) must be exact here. Indeed, we have tu^ = 
Tt 2 I/2mL 2 nJ AL = nu/uj i0) and 

(0) ALd^ 1 .7T 2 /. . 7T . 
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Thus the approximate relation (5.10) is valid here. It is also not difficult to 
check after a little algebra that in the exact relation (5.8) all higher powers 
of u cancel each other and the linear approximation is exact. 

If u = L is not extremely small, I ceases to be an adiabatic invariant and 
the system becomes nonintegrable. If L(t) is varied nonadiabatically and 
periodically around Lq, the system (termed the Fermi model [7B]) can show 
steady energy growth and chaotic motion. 

The generalization of (5.10) to multidimensional integrable systems with 
actions I\,...,If is evident, as is the generalization to multiparameter systems, 
with parameters Ai, A2, .... 

The 0(A) corrections in have been derived by other methods [77] . 
but in different forms from (5.10). If we assume A(t) is periodic, with a long 
period T\, the 0(A) terms in when integrated over (0,T\), define the 

Hannay angle A9h, the classical analogue of the quantum Berry phase [77]. 
The Hannay angle is the extra (geometric) phase undergone by the system 
in one cycle of A(i), in addition to the more obvious dynamical phase do ob- 
tained by integrating the instantaneous frequency u(\(t)) = dE(I, \{t))/dl. 
In general one requires more than one A-type parameter to obtain a non- 
vanishing result for A9h , but Hannay and Berry and coworkers give a few 
examples [77J [THJ EHj where one parameter suffices. We discuss two examples 
with nonvanishing A##. 

5.2.1 Bead Sliding on Slowly Rotating Horizontal Ring 

Consider a bead of mass m sliding without friction on a closed planar wire 
loop which is slowly rotating with angular velocity Q in the horizontal plane. 
For a loop of arbitrary shape, Hannay and Berry [77] show that the Hannay 
angle A9h defined above is given by A6h = —8tt 2 A/L 2 , where A is the area 
and L the perimeter of the loop. For a circular loop rotating about an axis 
perpendicular to the plane of the loop and through the center, A8h is simply 
— 27r, as is clear intuitively: the "starting line" for the particle moves ahead 
by 27r while the particle rapidly slides around the loop. 

We consider a circular loop of radius a (Fig. 8), but with the vertical axis 
of rotation through a point O on the wire instead of the center C of the 
loop, to make the answer less obvious. The adiabatic parameter here is the 
prescribed orientation of the loop A(t), where A(t) = Q(t) is small. We do 
not at this stage assume Q(t) is constant. 

Using the horizontal x, y lab-frame (Fig.8), we can easily write the La- 
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Figure 8: Circular wire loop rotating slowly in the horizontal x, y plane with 
angular velocity Q = A about axis through O perpendicular to the x, y plane. 
At t — 0, the particle starts with 9 = and the rotating x' axis ( a diameter) 
coincides with the fixed axis x. 

grangian, which is just the lab kinetic energy of the particle. Using 9 in Fig. 8 
as the generalized coordinate we find 



T 1 2 

L = -ma 
2 



Q 2 + {Q + 9) 2 + 2Q{Q + 9)cos9\ . (5.11) 

From (5.11) we obtain the equation of motion 

8 + (1(1 + cos9) +tt 2 sin9 = . (5.12) 

If the loop rotates uniformly (Cl — 0), (5.12) is just the pendulum equation, 
as one expects intuitively. 

Because we are using a rotating coordinate system (6 is measured from 
x') L is not equal to K — V = K = \ma 2 Q 2 in this system, so that the 
Hamiltonian H is not equal to the energy (= ^ma 2 9 2 ) in this system either. 
To find H, we first calculate the canonical momentum p e , where 

QT 

p e = = ma 2 (tt + 6 + Qcos 6) , (5.13) 
89 

and then H = pg9 — L is 

H = ^ma 2 9 2 -ma 2 Q 2 (1 + cos 9) . (5.14) 
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In terms of the canonical variables 8, pe the Hamiltonian H is 



2 

H= -E0- P9 Q(1 + cos 8) + 0(n 2 ) , (5.15) 

where we drop 0(Q 2 ) terms since we assume Q <C 8. 

We cannot apply (5.10) as it stands, since H = E was used in the deriva- 
tion. The generalization to include cases where H ^ E is 

dH . • oH . d . 1 , 

where H is the average of H over the fast motion (here 8). Alternatively, we 
can use (5.10), with now = dH(I, X(0))/dI. From (5.15) we have 

2 

ff=-^-pA (5.17) 
2ma z 

since 8 varies rapidly and pg is nearly constant. 

We now assume Q is constant. For Q — 0, H is a function of the good 
action I = pe only (no 8 dependence), and when fl ^ we see pe remains a 
good action. Note that H is independent of angle A, so that the second term 
in (5.16) will not contribute to uj{j>q) = dH/dpe, which is 

The total phase undergone after a long return time = 2tt/Q is thus 

8 = (^- 2 )T n -2n . (5.19) 
ma z 

The first term is the dynamical phase, the second is the Hannay (geometric) 
phase. 



5.2.2 The Foucault Pendulum 

For a Foucault pendulum the Hannay angle is the angle by which the plane 
of the pendulum vibration shifts in one day, i.e. A0# = — 2n cos 8, where 8 is 
the co-latitude of the point on earth where the pendulum is located. This is a 
clockwise rotation in the northern hemisphere. This result for has been 
derived by Berry [75] and others jHO] by various arguments. There are also 
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some pre-Hannay purely geometric arguments for the Foucault pendulum 
precession Here we show the result follows simply from (5.16). 

We consider a pendulum of length L and bob mass m oscillating at a 
point on the earth's surface with co-latitude 9. We use axes x,y,z attached 
to the earth, with x,y in the horizontal plane and z vertical. When at rest in 
the x,y,z system, the pendulum hangs along the z-axis, and when in motion 
the bob oscillates very nearly in the x,y plane. Neglecting the z motion and 
assuming Q = Qe cos 8 is small, where Qe is the earth's angular velocity, we 
easily find [H2| the approximate horizontal equations of motion 

x + u^x- 2Qy = , y + uJoV + 2Vtx = , (5.20) 

where ujq = (g/L) l l 2 . We have kept the 0(Cl) Coriolis terms, but dropped 
the 0(Q 2 ) centrifugal terms, since Q is small compared to the pendulum 
frequency u . Here Q is the local vertical component of the earth's angular 
velocity. We are interested in particular in the orbit precession rate predicted 
by (5.20). 

It is easily checked that (5.20) follow from the Lagrangian 

L = -m(x 2 + y 2 ) — -mujQ^x 2 + y 2 ) + mQ(xy — yx) , (5-21) 

and of course (5.21) can be derived directly 

Changing to polar coordinates p, <p y i a x = p cos <ft, y = p sin <p changes L 
to the form 

L = -mp 2 + -mp 2 (fi 2 — -mu^p 2 + mQp 2 (p . (5.22) 

Two constants of the motion are obvious from (5.22): since L does not depend 
on 0, is a constant, where 

dL 

Pcf, = — - = mp 2 (j) + mp 2 Q ; (5.23) 

dcf) 

and because L does not depend on t explicitly, the Hamiltonian H = p p p + 
p^tp — L is a constant, where p p = dL/dp and 

H = ^mp 2 + ^mp 2 <fi 2 + ^muj 2 p 2 . (5.24) 

Two constants of the motion with two degrees of freedom imply that the 
system is integrable and that two good actions exist. 
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To express H in terms of the canonical variables p, 0, p p , p^ we use (5.23) 
and 

dL 

p p = — = mp (5.25) 



to get 

H = i + ^ + ~ np *' (5 ' 26) 



where again we drop an 0(Q 2 ) term. We wish to use (5.16) (not (5.10) here 
as H 7^ K + V) to find the frequencies u p = u(I p ,X) and = w(J^, A), 
where A is the orientational angle of the earth multiplied by cos 9. Note that 
for Q — A = 0, (5.26) is simply the Hamiltonian for an isotropic 2D harmonic 
oscillator. Calculation of the actions for the 2D oscillator, I p and 1$ = p^, 
is a standard exercise [El], and the Hamiltonian H(Q = 0) in terms of these 
actions is 

H(n = 0) = 2I p lu + \It\w Q . (5.27) 

The term Vlp^ = QI^ in (5.26) will not disturb the good actions, so that 
(5.26) becomes 

H = 2I p lu + \It\wQ - I^tt . (5.28) 
Application of (5.16) with I — I p and I = 1^ then gives 

ujp = 2u> , ojtf) = ±uo — Q , (5.29) 

where a positive (+) frequency in corresponds to a counterclockwise ro- 
tation. Note that the Foucault precession frequency u>f = — is true for 
orbits of any shape, e.g. near linear (the usual case), elliptical, or circular. 
The factor of 2 in u p in (5.29) is clear intuitively: for Q = 0, in one circuit 
of an elliptical orbit, p goes through two cycles while </> goes through one. 
Computing the phase of the pendulum 4>(Te) after one complete rotation of 
the earth (period Te) then gives 

(f)(T E ) = ±uj T e - fLT E = ±uJqT e - 2vr cos 9 . (5.30) 

The term ±uqTe is the dynamical phase (f)o and — 27rcos# is the geometric 
Hannay angle A0#. 

The examples in this section and the previous one led to a nonvanishing 
Hannay angle A9h- It is easy to show that the earlier example (particle in 
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a ID box with a slowly moving wall) gives a vanishing result. The same is 
true for a particle in a slowly rotating ID box (fixed walls). Golin [HS] has 
proved that, for ID systems with one adiabatic parameter A, one always has 
A9h = unless A undergoes a rotational excursion through 2tt. However the 
proof does not cover our rotating ID box example (which has a hard wall 
potential), since Golin (SHI assumes a smooth potential. The case of a slowly 
rotating 2D box of elliptical shape has been worked out by Koiller jHO], and 
gives a nonvanishing A6h depending in a complicated way on the area of the 
box. 

5.2.3 Nonintegrable Systems 

Nonintegrable systems have fewer than / independent constants of the mo- 
tion or good actions, where / is the number of degrees of freedom. For such 
systems, at least some of the possible motions are chaotic, and for these at 
least some of the good action variables do not exist. Hence we have not been 
able to proceed rigorously beyond the general relation (5.4) for chaotic mo- 
tions. (Despite the nonexistence of good actions / for nonintegrable systems, 
formal use of them and their adiabatic invariance have been used success- 
fully in semiclassical quantization schemes; see sections 2.3, 3.1, 4.1, 4.3 and 
ref. |T3] for a method based on the Maupertuis Principles for conservative 
systems, and see refs. jnUE21 for a method based on adiabatic switching.) 

As stated, we do not have rigorous arguments to draw consequences from 
the UMP (5.4) for nonintegrable systems, but some heuristic arguments can 
be given. For example, by choosing H(\(t)) to vary slowly and linearly with 
t, e.g. H(X(t)) = \(t)H (q,p) = (1 + \t)H (q,p), where A is small, we can 
derive the adiabatic invariance of the action / for a ID periodic system (ear- 
lier we assumed that I is an adiabatic invariant for such systems). Further, 
by considering long trajectories of multidimensional systems, which start and 
end in some small region in phase space, and which are not quite periodic 
in general, we can formally show that, in a Poincare cycle time or in the 
limit T — > oo the long-path action W is an adiabatic invariant for quasiperi- 
odic and chaotic trajectories of multidimensional systems. Ehrenfest jHSj has 
previously shown that W is an adiabatic invariant for multidimensional pe- 
riodic systems. For ergodic systems, the phase volume Q(E) is an adiabatic 
invariant (see next paragraph), and since the adiabatic invariant is appar- 
ently unique for such systems (Kasuga [HI]]), it should be possible to derive 
the adiabatic invariance of W from that of Q for these systems, but we have 
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been so far unable to do this [EH]- We note that counterparts of the Hannay 
angle have been established for adiabatically evolving nonintegrable [EE] and 
ergodic jHH] systems. 

For one limiting case, complete chaos (ergodic motion) where the energy 
is the only constant of the motion, and long trajectories typically cover es- 
sentially the whole energy surface, a rigorous adiabatic invariant has been 
found jnni EH E3] , i.e. the phase volume fl(E), enclosed by the energy surface 
H(q,p) = E: 

Q(E)=J dqdp&[E - H(q,p)], (5.31) 

where G is the step function, which is unity for positive argument. (For 
multidimensional systems, dqdp stands for dqidq 2 . . .dqfdpidp 2 . ■ -dp f .) 

If H(q,p, A(i)) depends on a slowly varying parameter X(t) so that E = 
E(X(t)), then Q(E(X(t))) = Sl(E(\(0))). The adiabatic invariance of Q(E) 
is intuitively plausible First consider the initial (t = 0) energy surface 
H(q,p, A(0)) = E(X(0)), and its image under an arbitrary time development 
(i.e. not necessarily slow). The final (image) surface at the time T is not 
in general a constant energy surface, and hence not in general a dynamical 
surface. However, by Liouville's theorem the volume enclosed by the final sur- 
face is equal to that enclosed by the initial surface. If the time-development 
from t = 0toi = Tis sufficiently slow then at each stage the surface which 
has evolved from the t — surface is a dynamical surface, by definition of 
adiabatic evolution, where allowed motions transform into allowed motions 
PH] . In this case the final t = T surface is a true constant energy dynamical 
surface, and hence Q(E(X(T))) = Q(E(X(0))) from Liouville's theorem . 

For ergodic systems, there are no good quantum numbers except for the 
energy-ordering quantum number N = 1,2,..., where E\ < E 2 < .... A 
rough but useful semiclassical quantization scheme for ergodic systems based 
on Q(E) is the Weyl rule ^ 

Q(E N ) = Nh f , (5.32) 

which implicitly defines En- This is just a restatement of the early rule 
due to Planck j^H] that each quantum state occupies a volume of size in 
classical phase space. In practice, various improvements such as replacing 
N by (N + a) as in (2.4) are used jHS]- More subtle improvements can be 
obtained from Gutzwiller's trace formula [2E] • These relate fluctuations from 
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the smooth dependence of E on N predicted by (5.32) to the classical pe- 
riodic orbits. In cases such as the x 2 y 2 potential, where the classical Q(E) 
diverges [§E| due to the infinite channels of the potential (see Fig.l), one 
can estimate the smooth part of N(E) ( the cumulative number of states up 
to E) quantum mechanically jHZI at low energy by an asymptotic expansion 
JHEIj with leading terms e 3//2 loge and £ 3//2 , where e = E(m 2 /h C) 1 / 3 is di- 
mensionless. Application of the generalization of (5.32), i.e. N(E^) = N, 
then gives reasonable estimates of the low-lying state energies [3*61 198j. 



5.3 The Classical Hellmann - Feynman Theorem 

Another consequence of the UMP (5.4) is the classical version of the Hellmann 
- Feynman theorem 

dE(I, A) dH(q,p,X) ^ 

~^x~ ~ K — d\ — >I ■ (5 - 33) 

Here H(q,p, A) depends on a parameter A, so that the energy E(I, A) depends 
on A, as well as on the (constant) action J, and < ... >j denotes a time- 
average over a period T(I) at fixed action /. We assume a one-dimensional 
periodic system for simplicity. 

The quantum version of the theorem |lUUj is nowadays much better 
known, i.e. 

9E n (\) dH(X) 

dX =< — >»' (5 ' 34) 

where < ... > n =< n\...\n >, and where the Hamiltonian operator H(X) 
depends on some parameter A, so that the energy eigenvalue E n (X) and 
eigenfunction ip n (X) do also, because of 

H(X)M^=E n (X)M^)- (5-35) 

Proof of (5.34) is simple; we apply d/dX to both sides of the relation E n (X) = 
(ip n {X)\H (X)\ip n (X)) , and note that the two terms involving dip n (X)/dX add to 
zero because of (5.35), the hermiticity of H(X), and the fact that we assume 
(V„(A)JVn(A)) = 1 for all A. 

It is clear that (5.33) is the classical limit of (5.34), since as h — > (or 
n — > oo) (...)„ = (n\...\n) approaches a time average over a period, and fixing 
n corresponds to fixing /, since I ~ nh. To derive (5.33) from (5.4) we 
consider H(X(t)) to evolve adiabatically, so that A and T in (5.4) can be 
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related. The derivation is intricate, however, and it is simpler to proceed 
instead from the RMP (2.3) for conservative systems. We apply (2.3) to two 
periodic trajectories of H(\), the true one with period T(A), and a virtual 
one taken as the true trajectory of if (A + AA) with period T(A + AA) . Since 
first-order variations in E =< H(X) > vanish (at fixed W = 2ttI), we have 

< H(X) > A+A A=< H(X) > A +0(AA 2 ), (5.36) 

where < ... >a is an average over the trajectory with period T(A), and it is 
understood that I is the same for both < ... >a and < ... >a+aa- From (5.36) 
we get for AE =< H(X + AA) >a+aa - < H(X) >a to O(AA) accuracy 

AE =< H(X + AA) >a+aa - < H(X) > a+ aa=< H(X + AA) - H(X) > a+ aa • 

Again neglecting 0(AA 2 ) terms we have 

AE =< H(X + AA) — H(X) >\ . (5.37) 

Eq(5.37) yields dE/dX =< dH(X)/dX > A , and since E = E for the true 
trajectory we get dE/dX =< dH{X)/dX >\, the Hellmann-Feynman relation 
(5.33) apart from notation. 

As a simple example of (5.33) consider a ID harmonic oscillator, with 
Hamiltonian H(x,p) = p 2 /2m + kx 2 /2 and E(I) = Itu, where to = (k/m) 1 ^ 2 
is the angular frequency. Choose A = k, the force constant. The left side 
of (5.33) is E/2k and the right side is < x 2 /2 >=< V > jk. Thus we have 
< V >= E/2, as is well known from the virial theorem [Sj. If we choose 
A = m, we find similarly < K >= E/2. 

As a second example, consider the bound motion of a ID system with 
Hamiltonian H(x,p) = p 2 /2m + V(x). We change variables x,p — > x',p' via 
the transformation x' = Ax, p' = p/X, where A is an arbitrary scale factor. 
This transformation is canonical and preserves the Hamiltonian |lUlj . i.e. 
H(x,p) = H\(x',p f ), where 

H x (x',p') = X 2 ^ + V(^) . (5.38) 
2m X 

The energy of a given state of motion is unchanged by the change from H to 
H\ for any value of A. We now use (5.38) in the Hellmann-Feynman relation 
(5.33). The left hand side vanishes. The right hand side derivative is 

dH x Xp' 2 x> y , 2 R x y , 
dX m X 2 X X 
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where V = dV(x)/dx, and we have reverted to the original variables in the 
second form. Hence we get 



< 2K >=< xY > , 



(5.39) 



i.e. the virial theorem for an arbitrary potential V(x). 

6 Variational Principles in Relativistic Me- 
chanics 

Hamilton's Principle is widely used in Classical Relativistic Mechanics and 
in the Classical Theory of Electromagnetic and Gravitational Fields to derive 
covariant equations of motion (see, e.g. |ll)2p . As for Maupertuis' Principle, 
it is widely believed that it is not well suited for that purpose because when 
operating with the energy one loses explicit covariance. In this section we 
demonstrate, first, how to use Maupertuis' Principle for the derivation of the 
covariant equations of motion and, the second, how to use it for the solution 
of specific problems. 

Consider a relativistic particle with mass m and charge e moving in an 
arbitrary external electromagnetic field, described by a four-potential with 
contravariant components A a = (Aq,Ai) = ((f), Ai), and covariant compo- 
nents A a = (A ,—Ai) = ((p,— Ai), where and Ai (for i = 1,2,3) are the 
usual scalar and vector potentials. Hamilton's Action for this system can be 
written in the Lorentz invariant form |1()2[ ll()3j 



where we use the sign of Lanczos |lU3j for S, opposite to that of Landau 
and Lifshitz |lU2j . The sign can be chosen arbitrarily; our choice allows us 
to use below the orthodox definitions p = dL/dv and H = pv — L. The 
four- dimensional path runs from the initial point xa = (x% x\, x\, x A ) to 
the final point xb = (x B ,x B ,x B ,x B ) in four- dimensional space-time with 
corresponding proper times sa and sb- Here ds is the infinitesimal interval 
of the path (or of the proper time) ds 2 = dx a dx a = g a pdx a dx^ = dx\ — dx\ = 
dt 2 — dxj, the metric has signature (+, — , — , — ) and we use the summation 
convention and units with c = 1, where c is the velocity of light. S is itself 
not gauge invariant, but a gauge transformation A a — > A a + df /dx a ( for 




(6.1) 
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arbitrary /) adds only constant boundary points terms to S, so that 5S is 
unchanged. The HP is thus gauge invariant. 

If we introduce a parameter r along the four-dimensional path (proper 
time s along the true or any particular virtual path are valid choices) the 
action S can be rewritten in the form 

S = mJ [v a v a ] l/2 dr + e J A a v a dr, (6.2) 

where v a = dx a /dr. In general, r is both frame-independent and path- 
independent. A path- independent parameter is here invaluable for variational 
purposes, since we then do not have to vary the parameter when we vary the 
path. After the variations have been performed, we can then choose t = s, 
etc. The limits of the integrals are the invariants ta and t#. With respect 
to parameter r we get the covariant Lagrangian 

L = m[v a v a } 1/2 + ev a A a , (6.3) 
and corresponding conjugate momenta 

„ dL 



P = 7T~ 

and Hamiltonian 



mv a /[v0V fi ] 1/2 + eA a , (6.4) 



H = p a v a - L = 0. (6.5) 

Thus in this particular covariant treatment the Hamiltonian is trivial. 
The dynamics is hidden in a constraint. One can see that the four conjugate 
momenta are not independent variables but satisfy the constraint 

( Pa -eA a )(p a -eA a ) = m 2 , (6.6) 

which is obvious from (6.4). This constraint can be taken into account with 
the help of the Lagrange multiplier method, i.e. with the help of an effective 
Hamiltonian [103 EH 

H eff = X[(p a - eA a )(p a - eA a ) - m 2 ]. (6.7) 

We easily find that the Hamilton equations of motion with H e ff are equiva- 
lent to 

A = l/2m 
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(for t = s, the true trajectory proper time ) and to the covariant Lorentz 
equations of motion for a charged particle in an external field: 

dl> a a 

m— = eF af3 v p , 

_ dA p 8A a 
Fa '-d^~ drf' (6 - 8) 
i.e. the same equations found from the Lagrangian (6.3) (see, e.g. 

BED- 

On the other hand the Hamilton equations with H e ff are equivalent to 
the MP or to the RMP with the same H e ff . Thus one can use any of the 
four Variational Principles (HP, RHP, GMP, RMP) either with Lagrangian 
(6.3) or with effective Hamiltonian (6.7) for the derivation of the covariant 
equations of motion. It is to be noted that in the covariant formulation of 
relativistic mechanics, the Hamiltonian (a Lorentz scalar) is not equal to the 
energy ( the time component of a four- vector) . Thus E in the Maupertuis 
Principle must be interpreted as 

1 TB 

E = H eff = —j drH eff , (6.9) 

T A 

where At = r E — T4. By contrast, in the noncovariant formulation discussed 
below, the Hamiltonian is equal to the energy K + V + m . 

There are other, related, formal difficulties arising from (6.5). This con- 
dition (H = 0) is due to the fact that the Lagrangian L in (6.3) is homoge- 
neous and of first degree in v a (i.e. essentially linear in velocities v a ), so that 
v a dL/dv a = L and hence the difference v a dL/dv a — L (the Hamiltonian) 
vanishes |lU4j . Since H = 0, the Legendre transform relation (6.5) between 
L and H simplifies to 

L = p a v a } (6.10) 

and hence the Legendre transform relation between W and S, i.e. W — S = 
HAt (where H = E), simplifies to 

TB %B 

S = W = J Pa v a dr = J Pa dx a . (6.11) 

Because H = 0, direct application of the GMP (5W) E=0 = would be cum- 
bersome, and direct application of the RMP (SE) W = would be impossible. 



50 



All these difficulties can be traced to the well studied difficulty HU3| 
105, lOGJ ll(J7j of applying an action principle with kinematic |l(J9j non- 
holonomic constraints, such as (6.6) for the momenta, or for the velocities 
v a v a = 1 for proper time s, or v a v a = (ds/dr) 2 for arbitrary parameter r. 

One way out of the difficulties is to introduce the constraint with a La- 
grange multiplier for either the Lagrangian or Hamiltonian (as is done above 
for H e ff). Another way out is to exploit the nonuniqueness of the Lagrangian. 
For example |lU5j . the covariant Lagrangian L (6.3) is equivalent to the co- 
variant Lagrangian L , where (with r as parameter) 

r j I 

U = \{r)—v a v a + ev a A a , (6.12) 

where A(r) = dr/ds. It is easily verified that (6.12) generates the same 
Euler-Lagrange equations of motion as (6.3). For r = s (the proper time of 
the path in question), (6.12) simplifies to 

L' = -v a v a + ev a A a } (6.13) 

which is often employed Q-06J. The advantage of (6.12) and (6.13) is 
that they are nonhomogeneous in v a , and hence will generate a nontrivial 
corresponding Hamiltonian H' . For (6.13), for example, we find 

H' = ^-(p a - eA a )(p a - eA a ). (6.14) 
2m 

Note that H' (6.14) and H e ff (6.7) (for r = s) differ only by a constant. Use 
of the covariant V and H' in place of L and H removes the difficulties of 
using all four variational principles (HP, RHP, GMP, RMP). For L and H, 
only the first two can be used. Note again that E is computed from H' as in 
(6.9), with H' in place of H e ff. 

Now we come to the question of the solution of specific physical problems. 
In this case explicit covariance usually is not of great importance and it is 
more useful to find an appropriate convenient frame. For any given frame 
we can choose r — t (the coordinate time) and obtain the complete set of 
nontrivial quantities 

T 

L = -m(l - v 2 ) 1/2 - e(f) + ev • A, S = J Ldt, 

o 
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raw 

P = (1-^2)1/2 + eA ' 



T T 2 



p-vdt = J dt[ + ev • A], (6.15) 

^ ' 



The expression for L in (6.15) is derived from (6.3) or (6.12) by choosing 
r = t, except that we have used the freedom to change the sign of L to agree 
with the standard choice in the noncovariant case 0. Again we have chosen 
£a = and ts — T. Note that the nonrelativistic form W = J T 2Kdt for the 
Maupertius action is not valid here, and also note that the Hamiltonian H 
is equal to the energy K + ecp + m, despite the fact that L ^ K — V. 

As a result the complete set of four variational principles is available 
in any particular frame. It is matter of taste or practical convenience to 
use one principle or another. From a practical point of view there is no 
difference in the usage of variational principles in nonrelativistic mechanics 
and in relativistic mechanics in a given frame. 

As an example of the usage of the GMP for relativistic mechanics consider 
the motion of an electron in a uniform magnetic field B = (0,0,5). A 
corresponding vector potential in the Coulomb gauge (V ■ A = 0) is A = 
B x r/2, or A = (-yB/2,xB/2,0) [HE]. The solution of (6.8) in this case 
is well known - the electron moves with constant velocity along the direction 
of the field (along the z-axis in our case) and along a circular orbit in the xy 
plane with the cyclotron frequency 

(6.16) 

where E is the total energy of the particle including the kinetic energy of the 
motion in the z -direction. 

To calculate this frequency from the GMP we take as a trial trajectory 

x = A cos uit, y = A smut, v z = 0, 

where u and A are free parameters. After some simple algebra we get 
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(6.17) 



In this case it seems that an economical method is to use the GMP (2.2). 
For this purpose one has to rewrite W as a function of E 

W^[l-{^f][uE-^Bl (6.18) 

and to calculate the variation SW at fixed E. The GMP (2.2) translates here 
to (dW/duj)^ = 0. One then easily finds the analytical solution 

\e\B 



The exact solution corresponds to the case when E = E. For our trial 
trajectory, we happen to have E = E and the approximate solution coincides 
with the exact one ( similar to the case of a harmonic oscillator). 

This simple example demonstrates how to use variational principles for 
relativistic systems. 



7 Classical Limit of Quantum Variational Prin- 
ciples 

In section 2.1 we mentioned that the classical limit of the Schrodinger Quan- 
tum Variational Principle is the Reciprocal Maupertuis Principle for periodic 
or other steady-state motions, i.e. 

(S^l-O^^l-O, CM) 

where H is the Hamiltonian operator corresponding to the classical Hamilto- 
nian H(q,p). Here T is the period for periodic motions, or T — > oo for other 
steady-state motions (quasiperiodic or chaotic). Solution of the left-hand 
side yields |118j the stationary state ip n with the energy eigenvalue E n , i.e. 

Hil> n = E n ij n , (7.2) 
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the time-independent Schrodinger equation. Details of the proof of (7.1) 
(which is straightforward only for integrable systems) are given in [T2*] . In 
ref . |l!9j it is explained how Schrodinger, in his first paper on wave mechanics 
in 1926, in essence worked backwards, i.e. starting from the right hand side of 
(7.1), he "derived" the left hand side, and thereby discovered wave mechanics. 
In his second paper on wave mechanics in 1926 he retracted this argument 
because his version of the right hand side of (7.1) was not quite correct. In 
this second paper he presented an alternative argument as a basis for wave 
mechanics (now in the standard textbooks) based on the analogy between 
geometric and wave optics on the one hand, and particle and wave mechanics 
on the other. 

We consider now the generalization of (7.1) to nonstationary states i/)(t), 
where the true dynamical states satisfy the time-dependent Schrodinger equa- 
tion 

ihdtif>(t) = Hip{t), (7.3) 

where dt = d/dt, and the initial state ip(0) is given. For notational simplicity 
we consider a single particle in three dimensions, so that ip = ip(r, t), and 

h 2 

H = K + V = -—V 2 + V(r), (7.4) 
2m 

where K = p 2 /2m and p = (h/i)V. The arguments which follow are also 
valid for nonconservative systems, where H = H(t). 

We shall consider here the quantum time-dependent variational principle 
due to Frenkel and Dirac |12(J| I121| I122j for the wave function i/j(t), but we 
note the existence of other quantum time-dependent variational principles, 
e.g. for the time-evolution operator and transition amplitudes [123J, for the 
Heisenberg operator version |124j and for the matrix mechanics version |125j 
of the Hamilton Principle, for the density matrix [126] , for the mean value of 
an arbitrary observable |127j . and for time-dependent density functional the- 
ory |128j . ( For completeness, we note that there are also other |13Uj quantum 
time-independent variational principles, in addition to that of Schrodinger.) 

The original Frenkel-Dirac variational principle is |12Uj 

< 8i)\H - i%d t \ijj >= 0, (7.5) 

for arbitrary variations 5ip(r, t) around the true dynamical state ^(r, t). Since 
(7.5) is true for all possible variations 5ip, then (H — ihd t )ip = for the true 
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states, as in (7.3). By adding (7.5) and its complex conjugate and using the 
hermiticity of H we can transform (7.5) to [129J 

5 < i)\H-ihd t \^ >= -ihd t < %j}\5il) > . (7.6) 

Integrating over an arbitrary but fixed time interval (0, T), and assuming 
if)(t — 0) is fixed, we find |131| from (7.6) 

S [ T < ip\H - %W t \i\) >dt = -ih < ipmmiT) > . (7.7) 

J 

Note that (7.7) is not in the traditional form for a variational principle 
S J Ldt = 0, because of the right hand side. It is more like the unconstrained 
principles discussed earlier. Some authors drop the right hand side without 
comment |132| . but it is required if arbitrary variations 5ip in ip are allowed; 
otherwise (7.7) will not generate the Schrodinger equation (7.3) and its com- 
plex conjugate jl29| . It is possible to obtain a vanishing right-hand side by 
restricting the variations 5ip in various ways |133j . e.g. requiring the trial 
wave functions if)(t) to be normalized at all times. Another way is to write 
the right hand side as —ih < 8ip(T)\ip(T) >*, formally regard ip and its 
complex conjugate ip* as independent (see below), and assume ip*(T) (as 
well as V'(O)) is fixed. However, we find it simplest to retain the general form 
(7.7), with no constraints on if)(t) except for the fixing of the initial condition 
ip{0). (When dealing explicitly with bound states later, we assume boundary 
conditions on ip(r,t).) 

By writing out the scalar products < A\B >= J drA(r)*B(r) in (7.7) and 
using the hermiticity of p = (h/i)'V, we can write (7.7) in field theory form 

5 £ dt J drC = -ih J driP(T)*5ip(T) , (7.8) 
where the Lagrangian density C is 

h 2 

£(^^ = ;HWf + W| 2 -^*<9^ . (7.9) 
2m 

Since ip is a complex function, whose real and imaginary parts can be varied 
independently, it is well known |134j that tp and ip* can be regarded as 
independent for variational purposes. These two variations in (7.8) with 
£ given by (7.9) will generate the Schrodinger equation and its complex 
conjugate. 
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Following Feynman |132| . we change variables —> A, S (or p,S) in 

the Lagrangian density C(if),if)*), where A(r,t) is the amplitude and S(r,t) 
the phase of the wave function, via the Madelung transformation 



i/j = Ae lS/n = ^pe iS/n , (7.10) 
where p = A 2 = \ip\ 2 is the probability density. In terms of p and S we have 

£(P ' S) = P 2m^ + 2^ {V ^ )2 + PV + p9tS ~ ~2 dtp - (7 - 11} 
In terms of p, S the boundary term in (7.8) is 

< ^(T)\8^(T) >= J dv{p^)5{p~{T)+ l -Jdvp{T)5S{T) . (7.12) 

We now take the classical (h — > 0) limit of (7.8). Using (7.11) and (7.12) 

we get 

5 [ T dt [ dr P V^- + V + d t S] = ( drp(T)5S(T) . (7.13) 
Jo J 2m J 

We have assumed that in the limit H — > 0, p and S approach well defined lim- 
its ( which we continue to denote by p and S). In particular we assume that 
Vy/p does not contain an 0(fT l ) term |136j . We vary p and S independently 
in (7.13). Varying p, with S fixed, gives 

J J dtdvbpl^^- + V + d t S] = 0. 
Since Sp(r,t) is arbitrary, the bracket [...] must vanish, so that 

d t S + ^^ + V = . (7.14) 
2m 

This is the classical Hamilton- Jacobi equation |138j for the classical action 
S(r,t) = f*Q Ldt' along true paths. 

Next we vary S in (7.13) with p fixed. We integrate the 5(VS) 2 term by 
parts using 

/ drp(V5S) -VS = - J dr5SV ■ (pVS) , (7.15) 
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where we drop a surface term by assuming p — > for r — > oo. This assumes 
the corresponding quantum state has finite norm. We integrate the S(pd t S) 
term by parts using 

[ T dtpd t 5S = p(T)SS(T) - [ T dtSSd t p , (7.16) 
Jo Jo 

where we use 8S(t = 0) = which follows from the fixing of ip(t = 0). Using 
(7.15) and (7.16) in (7.13) thus gives 

/ dt / rfr[V ■ (p— ) + d t p]5S = 0. 
Jo J m 

Since this is true for arbitrary SS(r,t), we get 

d t p + V-(p— )=0 . (7.17) 

m 

This is the continuity equation, reflecting conservation of total probability, 
since |138| V5" = p = mv is the momentum of the particle. Note that if the 
right hand side of (7.13) were missing, we would not get the correct continu- 
ity equation. In contrast to what we have done here, Feynman |132j keeps 
the H terms in (7.11), and the variational principle then leads to quantum 
generalizations |137j of the Hamilton - Jacobi (7.14) and continuity (7.17) 
equations. 

Thus we see that the classical limit of the quantum time-dependent vari- 
ational principle is the classical variational principle for the Hamilton- Jacobi 
and continuity equations. Interestingly the classical limit of the quantum 
time-dependent variational principle leads to a classical variational principle 
for the Hamilton action S along the true trajectories, whereas the classical 
limit of the quantum time-independent variational principle in the form of 
the left side of (7.1) leads to a classical variational principle for the trajecto- 
ries themselves (RMP -the right side of (7.1)). An alternative form |139j of 
the quantum time-independent variational principle, 

5 < tjj\H - E\i) >= 0, (7.18) 

will generate in the classical limit a variational principle for the Maupertuis 
(or time-independent) action W(r) along the true paths. In (7.18) E is a con- 
stant Lagrange multiplier introduced to remove the normalization constraint 
< >= 1. The extremizing ^(r)'s which satisfy (7.18) obey 

(#-£)V = 0, (7.19) 
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so that E is the energy eigenvalue corresponding to the eigenfunction ip. 

To derive the classical limit of (7.18) we proceed much as before. (We 
assume H is time- independent in this section.) We first apply the Madelung 
transformation to ip(r), i.e. 

V>(r) = A(r)e iW{r)/h , (7.20) 

where A(r) is the amplitude and W(r) the phase. Following similar steps as 
above, we find (7.18) becomes in the classical (h — > 0) limit 

5jdrp(r)[^^ + V-E} = , (7.21) 

where p(r) = A(r) 2 = |^(r)| 2 . 

Varying p(r) in (7.21) with W{v) fixed gives 

^E)1 + V -E = , (7.22) 
2m 

the time-independent Hamilton- Jacobi equation |138j for the Maupertuis ac- 
tion W(r) = / r r Q p • dr along true paths. If we vary W(r) in (7.21) with fixed 
p(r) we get the continuity equation 

vw 

V ■ p— = 0, (7.23) 

m 

since |138j VW = p = mv is the momentum of the particle. 

In ref. |T21, we show how the classical limit of (7.18) can be taken in a 
different way (using standard ID WKB wave functions) and for a ID periodic 
orbit leads to 5(S + ET) = 0, i.e. the UHP for a periodic orbit of period T. 

We note that, just as for the right hand side of (7.1), the derivation of the 
classical limits (7.14) and (7.22) is somewhat formal for the chaotic motions 
of nonintegrable systems, since classical values of S and W may not exist 
as well defined functions for such systems. Unlike the nonrelativistic limit 
(v/c — > 0), which is regular, the limit h/S — > or Ti/W — ^ is singular, 
which makes the classical limit subtle and interesting |14(Jj . 

8 Summary and Conclusions 

We summarize the two groups of classical variational principles reviewed in 
this paper. First we have the Maupertuis principles: 

6W-T5E = 0, (8W)e = 0, (5E) w = 0. (8.1) 
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The first of these is the unconstrained form (UMP), the second the general 
Maupertuis principle (GMP), and the third is the reciprocal Maupertuis 
principle (RMP). By Legendre transformation of the principles in (8.1) we 
obtain the Hamilton principles: 

6S + E6T = 0, (8S) T = 0, (8T) s = 0. (8.2) 

The first is the unconstrained form (UHP), the second the Hamilton principle 
(HP), and the third the reciprocal Hamilton principle (RHP). In all these 
principles, the end-positions of the trajectory q& and qs are held fixed, and 
the notation indicates what else is fixed and what is varied. 

We have given a number of examples showing how to use these principles 
to solve practical problems, focussing particularly on the three Maupertuis 
principles. Both classical and semiclassical applications have been given, and 
some comparisons are made with the results from quantum variational prin- 
ciple calculations. The classical limits of two quantum variational principles 
(time-dependent and time-independent) are given. In particular, the RMP is 
the classical limit of the Schrodinger time-independent quantum variational 
principle. The RMP is thereby naturally adapted to semiclassical applica- 
tions 

In the applications, the variational principles ( both classical and quan- 
tum) are solved by the direct variational method, where one guesses a trial 
solution containing adjustable parameters, and adjusts the parameters to 
satisfy the variational principle. The Euler-Lagrange differential equation is 
not required in this method. When we choose a simple trial solution, the 
results can be found simply and analytically. 

Using the direct variational method the RMP and UMP have been used 
to estimate semiclassically the energy levels in several examples. The lev- 
els E nitn2 , ... can be expressed in analytical form, and the accuracy is quite 
acceptable in most cases, even for nonintegrable systems which are usually 
treated by much more complicated semiclassical methods, e.g. the Gutzwiller 
trace formula j2H]; or quantizing the Birkhoff-Gustavson normal form |141j . 
It would be useful to have comparisons of the results from the Maupertuis 
Principle-based method with those from the latter two and other methods 
for a number of specific problems. The methods should be compared for ac- 
curacy of results and computational complexity. Some limited comparisons 
of the Maupertuis Principle and normal form methods have been done for 
the x 2 y 2 oscillator of sec. 4.1 in refs. and limited comparisons for 
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this model of results from the Maupertuis Principle and the semiclassical 
adiabatic methods (35] are presented in Table 1 of sec. 4.1.1, but much more 
extensive comparisons would be of interest. 

Variational principles appeal to various individuals for differing reasons, 
both aesthetic |142j and practical. We have stressed the utility for classical 
and semiclassical particle dynamics calculations in this review. Engineers 
too make heavy use of variational principles in solving practical problems of 
classical continuum mechanics |143j . On rare occasions variational principles 
have even led to new laws of physics [144J. 
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10 Appendix I. Variational Notation and Re- 
ciprocal Principles 

"... since 5y is the variation of y with x kept constant, we could conveniently 
write it as (Sy) x , in accordance with a convention used in thermodynamics... 
It is curious that in spite of the obvious need in partial differentiation for 
precise statement of what is being kept constant, such statement is not em- 
bodied in customary notation of pure mathematics, though it is provided in 
thermodynamics, ..." 

H.Jeffreys and B.S.Jeffreys 
Methods of Mathematical Physics 
Cambridge U. P., 1946. 

Variational problems ask for the function which makes stationary an inte- 
gral involving the function and its derivatives. Usually there is a constraint 
which is fixed. For example, we are to find a function y(x), describing a 
closed planar curve, of fixed length L, which maximizes the area A enclosed 
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by the curve y(x). A and L are functional ofy(x): A = A[y(x)], L = L[y(x)]. 
For our purpose here, we do not need the precise form of A ot L. L is the 
constraint and A is the functional to be varied. We use the following notation 
to describe this variation 

(SA) L = 0, (All) 

where the constraint L is indicated as a subscript on the bracket enclosing 
the variation 5 A of the functional A. The variation 5 A is first-order in Sy, 
and is called the first variation when higher-order variations are considered 
15J . The constraint L is denoted explicitly in (Al.l) because it is the most 
important, or the one we wish to focus on. There may be other constraints 
which are left implicit, e.g. that the planar curve pass through a given point 
(xb,Vb)- in this paper we always leave fixed end-position constraints on 
curves (closed or open) implicit. 

The solution of the problem (Al.l) is the function y(x), which satisfies a 
differential equation, the Euler-Lagrange equation. The reciprocal variational 
problem has the constraint L and the varied functional A interchanged. In 
other words, it is described by the equation 

(SL) A = 0. (A1.2) 

The two equations, (Al.l) and (Al.2), have the same solution y(x). The 
reciprocal variation is searching for the function y(x) which makes stationary 
(here minimizes) the length L of a closed curve with fixed enclosed area A. 
The common solution for both problems is a circle, as is clear intuitively. 
One way to prove reciprocity is to use a trick, due to Lagrange, called the 
method of Lagrange multipliers. The problem (Al.l) is first restated in 
unconstrained form as 

5A-X5L = 0, (A1.3) 

where the variations of both A and L are unconstrained and A is an unde- 
termined multiplier. To obtain (Al.2) from (Al.3), we simply specify that 
we now fix A (i.e. set 5A = 0) and vary L, i.e. (5L)a = 0. This proves the 
reciprocity theorem. 

Reciprocity is clear intuitively in the isoperimetric problem (Al.l). An 
example of reciprocity from physics is the Gibbs conditions for equilibrium in 
thermodynamics |117j . The reciprocal conditions (5S)u = and (SU)s = 
are both conditions for thermal equilibrium, where S is the entropy and U the 
internal energy. These express the facts that the state of thermal equilibrium 
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is that state which maximizes the entropy for fixed internal energy, and also 
is that state which minimizes the internal energy for fixed entropy. (The 
unconstrained form |145j is SU — T5S = 0, where the Lagrange multiplier is 
the temperature T. Here the implicit constraints in all cases are fixed system 
volume V and mole numbers N a ) Reciprocal VP's in mechanics have only 
recently been studied (T^HH] - see Sec. 2. 

This is all we shall need from the calculus of variations; for more details 
one should consult a text-book |14b'j . 

11 Appendix II. Equivalence of the Mauper- 
tuis and General Maupertuis Principles 

In Sec. 2.2, for conservative systems, we discussed the equivalence of the 
Hamilton Principle (HP) 

(5S) T = 0, (A2.1) 
with the General Maupertuis Principle (GMP) 

(5W) B = 0. (A2.2) 

Here we show directly the equivalence of the GMP (A2.2) with the original 
Maupertuis Principle (MP) 

(SW) E = 0. (42.3) 

In all these principles the end-positions qA and qs are held fixed. Since 
(A2.1), (A2.2) and (A2.3) all generate the same (Euler-Lagrange) equations 
of motion, it is clear they must be equivalent, and it is interesting to show 
directly the equivalence of (A2.2) and (A2.3). We must show that relaxing 
the constraint of fixed E in (A2.3) to fixed E in (A2.2) does not generate 
spurious or unphysical solutions. 

The constraint of fixed E = H(q,p) in (A2.3) applies along and between 
the trial trajectories. To relax this constraint, we use the method of Lagrange 
multipliers |147j . Since H(q,p) = const, we have 5H = , and hence 

A [ T SHdt = 0, 
Jo 

where A is an arbitrary constant. The trajectories all start at t — at fixed 
point A and end at t — T (which may vary with the trial trajectory) at fixed 
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point B. We release the constraint of fixed E in (A2.3) and choose A such 
that 

SW = X [ T 5Hdt, (42.4) 
Jo 

where 5 indicates a variation from the true trajectory value and 5H ^ in 
(A2.4). This can be rewritten as 

SW = X[S [ T Hdt - H(T)5T], (42.5) 
Jo 

where H(T) = H(q(T),p(T)), and where we recall that the end-time T is 
not fixed in (A2.3) . Recalling the definition E = J T Hdt/T enables us to 
write (A2.5) as 

SW = X5(TE) - XE5T, (42.6) 

where we have used the fact that H = E = const on the true trajectory. 
Expanding (A2.6) gives 

5W = XE5T + XT5E - XE5T. (42. 7) 

The first and third terms on the right hand side of (A2.7) cancel, since 
E = E = const on the true trajectory, as assumed in the MP (A2.3).This 
means 

SW = XT5E. (42.8) 

To see that A = 1 in (A2.8), we apply (A2.8) to the special case of two true 
trajectories, with actions W and W + dW and energies E and E + dE. For 
this case (A2.8) reads dW/dE = XT. Comparison of this with the standard 
relation j3] dW/dE = T gives A = 1. Hence we have 

5W = T5E, (42.9) 

which is the unconstrained Maupertuis Principle (UMP) (2.10). Fixing E in 
(A2.9) then gives (A2.2), the GMP. 

Thus (A2.2) follows from (A2.3). The converse, the derivation of (A2.3) 
from (A2.2) is easy; we simply restrict the trial trajectories in (A2.2) to those 
of fixed energy E, since the equations of motion following from (A2.2) imply 
energy conservation. Hence (A2.2) and (A2.3) are equivalent. 

We note once more that conservation of energy is assumed in the original 
MP (A2.3), whereas it is a consequence of the general MP (A2.2). 
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